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Abstract. We give a proof of universality in the bulk for orthogonal (/3 = 1) 
and symplectic (/3 = 4) ensembles of random matrices in the scaling limit for a 
class of weights w{x) = e~^^^^ where V is a polynomial, V(x) = K2mx'^"^ + ■ ■ ■ , 
K2m > 0. For such weights the associated equilibrium measure is supported on 
a single interval. The precise statement of our results is given in Theorem im 
below. For a proof of universality in the bulk for unitary ensembles (/3 = 2), 
for the same class of weights, see 

Our starting point is Widom's representation IWI of the orthogonal and 
symplectic correlation kernels in terms of the kernel arising in the unitary 
case (/3 = 2) plus a correction term which is constructed out of derivatives 
and integrals of orthonormal polynomials (OP's) {pj}j>o with respect to the 
weight w{x). The calculations in IWI in turn depend on the earlier work 
of Tracy and Widom TW2 . It turns out (see )W| and also Theorems 12.11 
and Tn\ belowl that only the OP's in the range j = N + 0(1), Af ^ oo, 
contribute to the correction term. In controlling this correction term, and 
hence proving universality for [3 = 1 and 4, the uniform Plancherel— Rotach 
type asymptotics for the OP's found in X)KMVZ2 play an important role, 
but there are significant new analytical difficulties that must be overcome 
which are not present in the case /3 = 2. We note that we do not use skew 
orthogonal polynomials. 



1. Introduction 

In this paper we will be concerned with ensembles of matrices {M} with proba- 
bility distributions 

(1.1) pNf3iM)dM ^ -^e~''^f'^^UM, 

Zn,I3 

for p = 1, 2 and 4, the so-called Orthogonal, Unitary and Symplectic ensembles, 
respectively (see jHl)- For f3 — 1, 2, 4, the ensemble consists oi N x N real 
symmetric matrices, N x N Hermitian matrices, and 2N x 2N Hermitian self- 
dual matrices, respectively. In general the potential V/3(x) is a real- valued function 
growing sufficiently rapidly as \x\ — > oo, but we will restrict our attention henceforth 
to Vg's which are polynomials, 

(1.2) Vpix) = K2m,/3X'^"' + ■ ■ ■ , K2m,/3>0. 

In 11.1(1 , dM denotes Lebesgue measure on the algebraically independent entries of 
M , and Z]\[,i3 is a normalization constant. The above terminology for /? = 1, 2 and 
4 reflects the fact that 1(1. l|l is invariant under conjugation of M, M ^ UMU~^, 
by orthogonal, unitary and unitary-symplectic matrices U. It follows from ((1.1(1 
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that the distribution of the eigenvalues xi, • • • , sat of M is given (see [M]) by 



N 



(1.3) PN^j3{xi,- ■ ■ ,Xn) = Yi - ^kfti^f^'^^l^ 

l<j<k<N j=l 

where again Z^^p is a normahzation constant (partition function). Here 



(1.4) Wi3{x) 



^-2Vp(x)^ /3 = 4. 



(The factor 2 in wp^i reflects the fact that the eigenvalues of self-dual Hermitian 
matrices come in pairs.) Let {Pj}j>o be the normalized orthogonal polynomials 
(OP's) on R with respect to the weight w = W/3=2, and define 4>j = •pjW^I'^ . Note 
that (</)j,(/)fe) = (5jfc where (•, •) denotes the standard inner product in L^(M). 

For the unitary matrix ensembles (J5 = 2) an important role is played by the 
kernel 

(1-5) KN{x,y) = KN,2{x,y) = ^ (t>k{x) (l)k{y) ■ 

k=0 

In particular the probability density (|1.3|l . the ^-point correlation function Rn,i,2 
and also the gap probability £^2(0; J) that a set J contains no eigenvalues, can all 
be expressed in terms of K^, see e.g. For example 

(1-6) RnjM^i,--- ,xi) ^ det{KN{xj,Xk))i<j,k<i- 

The Universality Conjecture, in our situation, states that the limiting statistical 
behavior of the eigenvalues xi, - ■ ■ ,xjv distributed according to the law H1.3|) . in 
the appropriate scale as A'^ — > oo, should be independent of the weight Wfs, and 
should depend only on the invariance properties of Vn^P: /? = 1, 2 or 4, mentioned 
above. Universality has been considered extensively in the physics literature, see 
e.g. |E?71IB^IHWellHe7) . 

The kernel Kn{x, y) can also be expressed via the Christoffel-Darboux formula 

(1.7) KN{x,y)=bN-i , 

x~y 

where b^-i is a coefficient in the three-term recurrence relation for OP's (see (|2.1|) 
below). In view of the preceding remarks it follows that in the case /3 = 2, the 
study of the large N behavior of P/v,2, and in particular the proof of universality, 
reduces to the asymptotic analysis of b^-i and the OP's PN+j with j = or — 1. 
By a fundamental observation of Fokas, Its and Kitaev [FoIKi] the OP's solve a 
Riemann-Hilbert problem (RHP) of a type that is amenable to the steepest descent 
method introduced by Dcift and Zhou in DZ' and further developed in |DVZj . 
In [DKMVZll iDKMVZ2 the authors analyzed the asymptotics of OP's for very 
general classes of weights. In particular they proved the Universality Conjecture 
in the case /3 = 2 for weights w{x) = e"^'-^^ where V{x) is a polynomial as above, 
and also for w{x) — e"^^*^^) where V{x) is real analytic and V{x)/\og \x\ +oo, 
as |x| — > oo. The (bulk) scaling limit ^ oo is described in terms of the so-called 
sine kernel Koo{x — y) where 
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For example jDKMVZ'2l Theorem 1.4], for w{x) — e Vix) polynomial, and 

for any ^ = 2, 3, • • • and r, yi, in a compact set, one has as ^ cxd 

(1.9) 

w^^kw "^"^''''(^ ^dk^' ^^dk^) - 2/^-))i<.,/c<,. 

The scale x = y/Kf^(0, 0) is chosen so that the expected number of eigenvalues per 
unit 2/-interval is one: This scaling is standard in Random Matrix Theory. Indeed 
for any Borel set B C M, 

(1.10) / Rn,i=i,2{x) dx = E{ number of eigenvalues in B }. 

J B 

Thus by H1.6|l XAr(0,0) = -Rjv, 1,2(0) gives the density of the expected number of 
eigenvalues near zero. In other words, in the appropriate scale, the large N behavior 
of the eigenvalues is universal (i.e. independent of V). Pioneering mathematical 
work on the Universality Conjecture was done in PS^ and for the case of quartic 
two-interval potential V{x) = N{x'^ — te^), t > (sufficiently) large, in jBIj . We 
note again that all these results apply only in the case (3 = 2. 

In the case /? = 1 and 4 the situation is more complicated. In place of (|1.5|l one 
must use 2x2 matrix kernels (see e.g. [CT1ITW2] ') 
(1.11) 

SnA^^v) iSN,iD){x,y) 
KN^x.y) = ( I ' ^ even, 

{eSNs){x,y) ~ \sgTL{x - y) SN,iiy,x) 



and 

(1.12) KNAix,y) = ^ 



SnMx^v) iSN,4D)ix.,y) 



{eSNA)ix,y) SnaIv^x) 

Here SN,i3{x,y), /? = 1,4, are certain scalar kernels (see (I1.17f) . (|1.18|) below), D 
denotes the differentiation operator, and e is the operator with kernel e{x, y) = 
^ sgn(a; — yY- Such matrix kernels were ffist introduced by Dyson |Dy2| in the con- 
text of circular ensembles with a view to computing correlation functions. Dyson's 
approach was extended to Hermitian ensembles, ffist by Mehta fM2) for V{x) = x^, 
and then for more general weights by Mahoux and Mehta in |MaM| . A more direct 
and unifying approach to the results of Dyson-Mahoux-Mehta was given by Tracy 
and Widom in |T W2| . where formulae (|1.17(l . I|1.18|l below were derived. We see 
that once the kernels SN,ii{x^y) are known, then so are the other kernels in Kj^^fj. 
As in the case (3 = 2, the kernels i^Tv./S give rise to explicit formulae for RN,i.fj and 
Ep{Q; J). For example 

(1.13) RnaA^) = Ri,i3{x) = ^ti Kn,i3{x,x) 

and 

RN,2.^l3{x,y) = i {tiKN,p{x,x)){ti;:KN,p{y,y)) - ^tr {KN,j:j{x,y)KN,p{y, x)) , 

and so on, see |TW2| . As indicated above, formula (|l.lll) only applies to the case 
when N is even. When N is odd, there is a similar, but slightly more complicated, 
formula (see |AFNvM] ) . Throughout this paper, for /3 = 1, we will restrict our 

"'^We use the standard notation sgnx = 1, 0, —1 for x > 0, x = 0, x < 0, respectively. 
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attention to the case when N is even. We expect that the methods in this paper 
also extend to the case (3 — I, N odd, and we plan to consider this situation in a 
later publication. Of course, in situations where the asymptotics of has been 

analyzed (e.g. V{x) = cc^) for all N as N ^ oo, the limiting behavior of Rn,i,p=i 
is indeed seen to be independent of the parity of N (see e.g. |MllNWp . 

Let {qj {x)}j>o be any sequence of polynomials of exact degree j, qj (x) = qj.jx^ + 
• • • , qj,j 7^ 0. For j = 0, 1, 2, • • • , set 



(1.14) 



qj{x)wi{x), l3 = l 

q,{x){wi{x)Y'\ /3 = 4. 

Let Mjv,i denote the N x N matrix with entries 

(1.15) (Mwa),fc = (V'jM,eVfc4), 0<j,A:<A^-l, 
and let Mi\[^4 denote the 2N x 2N matrix with entries 

(1.16) (A/Ar,4),fe = (^j,4, i?V^fc,4), 0<j,k<2N-l, 

where again (•, •) denotes the standard real inner product on R. The matrices Mjv,i 
and Mn_4 are invertible (see e.g. |AvM[ (4.17), (4.20)]). Let fJ,N,i, fJ-NA denote 
the inverses of M^r^i, Mm, a respectively. With these notations wc have |TW2| the 
following formulae for Sm.p in (|l.ll(l . H1.12|l 

N-l 

(1-17) SN,i{x,y) = - ^ (/ijv,i)ife {eTpk,i){y) 

j,k=Q 
2N-1 

An essential feature of the above formulae is that the polynomials {qj} are arbitrary 
and we are free to choose them conveniently to facilitate the asymptotic analysis 
of Hl.ll|l . (|1.12l) as TV ^ oo. In view of the unitary case (3 = 2 described above, 
specific issues that must be addressed by the choice of the (?j's, are the following: 

(a) We need an analog of the Christoffel-Darboux formula to convert Sm,p into a 
form similar to (|1.7(l . which depends only on polynomials of high order as A'^ — *■ cxd. 

(b) By (a), the large A^ behavior of (|l.ll(l . H1.12|l becomes purely a question of 
the asymptotics of polynomials of high order. Thus we need to choose polynomials 
whose asymptotic behavior can be analyzed. 

(c) Finally the inverses of the matrices Mm,i and M^a must be controlled as 
A -> oo. 

Property (c) is new in the cases (3 = 1 and 4, and does not arise in the unitary 
case (3 = 2. Following |W|, we will choose the qj^s to be orthogonal polynomials 
with respect to an appropriate weight. It turns out that this choice addresses all 
three issues (a), (b), (c) simultaneously. 

Previously the study of the orthogonal and symplectic ensembles was often car- 
ried out via so-called skew orthogonal polynomials (SOP's) (see ^] for classical 
references and, e.g., [NWl lAFNvM. AvM. FNH for more recent work). Skew 
orthogonal polynomials are characterized by the property that (in our situation) 

Mm, (3 has the form of a direct sum of scalar multiples of the matrix ^ '^^ o) ' 
many classical weights wp, e.g., Wi3{x) = e^^ , the corresponding SOP's can be 
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computed explicitly jAFNvMl If] . The proof of universality, however, involves arbi- 
trary weights wp for which no such explicit formulae for the SOP's are available. As 
noted above, the crucial ingredient in the proof of universality for the (3 = 2 ensem- 
ble is the asymptotic analysis of the OP's corresponding to an arbitrary weight W2- 
For /3 = 1, 4, the problem of computing the asymptotics of SOP's raises challenging 
technical difficulties which have not yet been fully overcome (however see the work 

E). 

In order to state our main result we need more notation. For any m G N let 
V{x) be a polynomial of degree 2m 

(1.19) V{x) = K2mX^"' + ■■■ , K2^>0 

and let w{x) = wi3=2{x) = e~^(^^ as before. Let Pj{x), j > 0, denote the OP's with 
respect to w, and set 0j(a;) = Pj{x){w{x)y^^ , j > 0, as above. For /3 = 1,4 set 

(1.20) Vp{x) ^ ^V{x) 

and let be even. Then by H1.4|l . W4 = e^^^* = and wi = e~^^ = e"^/^. 
This ensures that for the choice qj = pj in H1.14|l 

(1.21) ipj.0=i{x) = ipj^0=4{x) = (t)j{x), 

which enables us in turn to handle Sn^i and 8^/2.4 in H1.17|l . (|1.18() simultaneously 
(see Remark 11.31 below). Henceforth and throughout the paper, Kn denotes the 
Christoffel-Darboux kernel H1.5() , H1.7|l constructed out of these functions 4>j . 

The main result we state here fTheorem ll.ll below') concerns universality for the 
kernels K^.p in the bulk scaling limit. Universality for other standard statistical 
quantities also holds, for example for the Z-point cluster functions and for the gap 
probability, see Corollaries 1 1 . 21 and 11.31 below . 

As in the case (3 = 2 above, we scale so that the expected number of eigenvalues 
per unit interval is one. Note that formula (|1.1U|I also holds for /3 = 1 and 4 and so 
Rn,i=i,i3{^) gives the density of the expected number of (simple) eigenvalues near 
zero. For /? = 1, 4, in view of ((T~T^ and ((T~TT|l . ifTT^ 

-Rjv,i.i(^) = SN,i{r, r) 

(1.22) 1 

RN/2,iA(r) = -j^SM/2A{r-,r). 

It follows from H4.154f) , H4.155|l below that as (even) N ^ 00 ioY r in a compact set 



SN,i{r,r)=KM{QM^ + 0{N-^'^)) 



(1.23) 
and hence 

i?jv.i,i(^) = ^w(0, 0)(1 + 0{N-^'^)) and 

(1-24) 1 

RN/2sAr) = -i^^(0,0)(l + O(A^-i/2)) 

give the scaling in Thcorcm ll . II below. 

Remark 1.1. The asymptotics in (|1.23|) . H1.24|l and in many related situations (see 
e.g. Remark il .21 below') remains true for r = o(A^^/'-^'"-'), but we do not need this 
fact in what follows. 
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It turns out that the off-diagonal elements in Kn.js scale differently as iV ^ oo. 
On the other hand, the statistics of the ensembles are clearly invariant under the 
conjugation 

for any scalar A — Xn,p- For example, this is obviously true for the cluster functions 

Tn I fs which have the form 

(1.25) 

where the sum is taken over ah permutations of {1, • • • , 1} (see |TW2I p. 816]), etc. 
Set 

Aat,! = {RNSA{r)f'^, Xna = (i?W/2,l,4(^))'^'- 

Theorem 1.1. Let j3 — 1 or A. For any V{x) of degree 2m as in (|1.19() define 
Vf3{x) and wp{x) as in (jl.2()|l . ()1.4|1 . Then for r, rj in a compact set, as (even) 
N oo 



(1.26) 



E^A ^ (r + ^,r+-f]- X(4)(e, ^) - 



whe 



(1.27) 




/o^-" (t) dt-^sgn{^-rj) 
K^{2i^-V)) ^(ifoo(2(C-r7))) 



''K^{2t)dt K^{2{7j-0) 
For /3 = 1, 4 and r, ^,r] in a compact set, we have as N ^ oo 



(1.28) 8n.,p = 



0(iV-i/4) o(7v-i/2) 



Remark 1.2. In the above theorem, we do not need to restrict r to a compact set. 
The theorem remains true for r = o(A^^/'^^™''). 

Remark 1.3. The choice in (|1.20l) implies that Vi{x) = Vi(x). Other authors 
make different choices. For example, in the Gaussian case |M], Mehta considers 
Vi{x) — 2Vi{x). Our choice insures that we can use the same auxiliary OP's in 
(|1.21|l for /3 = 1 and /3 = 4. On the other hand, Mehta's choice ensures that 
the density of states for Pm, pixi, - ■ ■ ,xn) in is the same for /3 = 1 and 

/3 = 4. Nevertheless we obtain precisely the same answer ()1.26|) in Theorem II .11 as 
computed in [M] for /? = 1 and /? = 4 in the Gaussian case. Of course, this is due 
to the fact that in the end we use the same scaling as in viz., the expected 
number of eigenvalues per unit interval is one. 
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From the formula for the cluster functions H1.25|l we immediately have the fol- 
lowing result. 

Corollary 1.2. Let (3 = 1 or 4. Let V be a polynomial of degree 2m and let K^^\ 
/? = 1, 4 be as in Theorem M.lX Then for [3=1 and Z = 2, 3, • • • we have, uniformly 
for r and , • ■ ■ , ^/ in a compact set 

um — 

Af^oo (A 

(1.29) 



lim — — -TN.j,i{r + — — ,---,r + 



For P = A, the same result is true provided we replace Tjv.z,! — > Tjv/2,i,4; Xn,i — 
AAr,4, and ^ 

Together with some additional estimates (see SectionOJ, Theorem ll . II also yields 
the following universality result for the gap probability. Recall that for a 2 x 2 
block operator A = {Aij)ij^i^2 with A\\,Aii in trace class and A\i,Ai\ Hilbert- 
Schmidt, the regularized 2-determinant (see e.g. [ST|) is defined by det2(/ + A) = 
det((/ + A)e-^)e*''(^"+^=^=^). 

Corollary 1.3. Let /5 = 1 or 4. Let V be a polynomial of degree 2m and let K^^\ 

(3 = 1,4, he as in Theorem M.lX Fix 6 > and r. Then 

(1.30) 



lim Prob { no eiqenvalues in \ r ^ — , r H — 1 ^' = \/ det'-'^V/ — K, 

A^-oo 1 V X%^^j XI M V 



(/3)n 



where Kg^^ denotes the operator with kernel K^^^{^,ri) acting in L^{—9,9). Here 
det''^-'(/ — Kg'^'^) is the regular determinant of the trace class operator Kg^\ but 
det''^-'(/ — Kg^^) is the regularized 2-determinant det2(/ — Kg^^) defined above. 

Remark 1.4. The regularized 2-determinant is needed for /? = 1 because the oper- 
ator with kernel ^ sgn(,^ — i]) is Hilbert-Schmidt but not trace class on L^{—9, 9) 
(see proof of Corollarv ll.3l in Section ISJ. 

We complete this introduction with a description of Widom's result jffij which 
is basic for our approach in this paper. Widom's method applies to general weights 
wp with the property that w'^/wp is a rational function. This property certainly 
holds for our weights as in (|1.4|l . (|1.2(l . Introduce the semi-infinite matrices 

(1-31) Doc = {{D(j)j,(l)k))j,k>0, ^ao = {{^4>j,4>k))jM>Q■ 

li follows from |TW1[ Section 6] that the matrix Dao is banded with bandwidth 
2n + 1 where 

(1.32) n = 2m-l. 

Thus (-Doo)j"fe = if \i — k\> n. Next, let N be greater than n, and introduce the 
following A'^-dependent n-column vectors 

$i(a;) = [<j)N-n{x), ■ ■ ■ ,cj)N-i{x))'^ 

(1.33) 

e$i(a;) = (e0Ar_„(x), • • • ,e(j}N-i{x)) 
e^2{x) = (e(/)Ar(a;), • • • , e^Ar+„_i(a;))^ 
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and the following 2n x 2n matrices consisting of four n x n blocks 

(Bii Bi2 
I ^ {{e(l)j,4'k))N-n<],k<N+r, 
B2I B22 

and 

(1.35) A = 



A12 \ / D12 
A21 I V -1321 



Dii D12 

where ( | = ((130^, (/)fe))iv-n<i,fc<Ar+n-i- Finally, set 

D21 D2: 

C = 





C12 ^ 




V C21 


C22 I 


H 




Note that 

C\\ = /„ + B12A21 = /„ — -Bi2-D2i- 
The main result in ^\ is the following pair of formulae for Sn^i and 8^/2^4 



(1.36) 
and 
(1.37) 



■{e<^,{yf,e<^2{yff 

SN/2Aix,y) = KN{x,y) + <i>2ixf -021 ■ e^i{y) 

+ <^>2ixf ■ D2iC^i^BnDi2 ■ e^2{y)- 



Observe that Sn_i and S'jv/2,4 £^re sums of the (3 = 2 kernel KN(x,y) together with 
correction terms that depend only on ipN+j for j G {~n^ ■ ■ ■ ,n — 1}. Together with 
the Christoffel-Darboux formula (|1.7II we see immediately that the representations 
(|1.36|l . H1.37|l address issue (a) above. Moreover, the ^Ar+j's are proportional to 
orthonormal polynomials PN+j > and hence can be evaluated asymptotically as iV — > 
00 by the methods in |DKMVZ2) . This addresses issue (b). Finally we see that 
issue (c), the question of the control of the inverses of the (large) matrices Mn,i, 
Mm A, reduces, via H1.36|l . H1.37|l to a question of controlling the inverses of matrices 
of fixed size 2n x 2n and n x n, respectively. It turns out that these fixed size 
matrices converge as iV ^ 00. The proof that in both cases the limiting matrices 
are invertible for any V of type 1)1.19(1 . is not obvious, and constitutes a significant 
part of the proof of Theorem II. II in this paper. 

In Section 12 we state a variety of auxiliary results which we use in the proof of 
Theorem 11.11 The results, some of which are of independent interest, are of two 
types: 

(i) asymptotics of {D(j)j\j+j, 4>N+k) and {e(j)N+j, 4>N+k) for fixed j, fc as — + 00; 

(ii) the equality of four specific determinants that arise in the analysis: This 
reduces the proof of the invertibility of the limiting matrices mentioned above to 
the proof of the invertibility of a single matrix Tm-i (see (|2.17|l below). 

In Sectional we prove Theorem ll.il fand hence Corollarv ll.2|l and Corollarv ll.3l 
assuming the validity of the results in Section |21( and In SectionQ] we prove the 
results of type (i) above. In Section we prove the equality of the determinants 
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in (ii). Finally, in Section El we prove that for any V{x) ~ K2mx'^™' + the 
determinant of Tm-i is nonzero. In the course of the proof in SectionlHl we need to 
estimate certain explicit integrals: A rigorous analysis of the error bounds in these 
estimates is given in an Appendix. 

Remark 1.5. After the second version of this paper was posted on the arXiv, 
A. Boutet de Monvel informed us of the work of A. Stojanovic f |Stll ISt2[ [StS]) 
in which universality (in the bulk and also at the edge) is proved for (3 = 1 and 
/3 = 4 in the even, quartic (two-interval) case considered previously by Bleher and 
Its |BI| for (3 = 2. The method uses a variant of the formulae in together with 
the asymptotics for OP's in jBI| . 

In ISt2| Remark 2.4], an interesting connection is noted between the problem 
of controlling limAr^oo det Cn (in our notation) and the problem of estimating 
partition functions Zjq^p (see (^31) as N ^ oo. Indeed from the formula Cn = In — 
B12D21 above and H2.18|l below, we see that det Cn = det e^r det Djy, where ejv, Dn 
are the leading NxN sections of the matrices £00 , Doo , respectively. But det e jv and 
detDpf can be expressed in terms of partition functions (see | AvM| ] and this leads 
directly to a formula for det Cn as a ratio of Z^^p's, /? = 1, 2, 4. Thus one may try to 
use the methods of statistical mechanics to control limjv^oo det Cn. The estimates 
in p] show that the partition functions Zj^^p have leading order asymptotics of 
the form e''"^ (i+°(i)) as iV ^ 00, and moreover, their combined contributions to 
detCii cancel to this order. Thus, in order to show that limAr^oo det Cn ^ 0, 
one needs higher order asymptotics for the Zjv,/3's. In one of our early attempts 
to prove universality, for (3 = 1 and 4 we indeed tried to derive such higher order 
asymptotics, but we were not successful. The statistical mechanical approach to 
prove limjv^co detCn 7^ 0, remains an unresolved, and intriguing, possibility. 

Notational remark: Throughout this paper c, C, C(m), ci, C2, • • • refer to con- 
stants independent of N, P. More specifically, the symbols c, C refer to generic 
constants, whose precise value may change from one inequality to another. The 
symbol cn however always refers to the iV-dependent constant 1)2. 3|l below. 
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2. Auxiliary results 

In order to evaluate the asymptotic behavior of the correction terms in (|1.36() 
and p. 37(1 . one first needs to determine the asymptotics of (I?(/)Ar+j, 0jv+fc) and 
{e(f>N+j,4'N+k) for fixed j,k e {— n, • • • ,n — 1}, as iV ^ 00. Such asymptotics 
are given by the next two theorems which are proved in Section 0] below. Let bj 
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be the coefficients in the three term recurrence relation |£z] satisfied by the OP's 
Pj corresponding to the weight w{x) — e^^'^^ (and hence also by the functions 

(2.1) xpj{x) ^ bj^ipj^i{x) + ajPj{x) + bjPj+i{x), j > 0, 

(&-i^0). 

Theorem 2.1. Let V(x) = K2mx'^™' + • ■ • ; '«2m > 0, and set n = 2m — 1 as before. 
Then the matrix Doo has 2n+ 1 bands and is asymptotic to the product of a diagonal 
matrix and a Toeplitz matrix. More precisely, for any fixed j,k G li, as N oo 
(2.2) 

{D(l)N+j,4'N+k) 

{0, j - A; = or |j - fc| > 71+ 1 

sgn(j - fc)((„_|j"fe|)/2) +o(l), |j-fc| = 1,3,-- - ,n 

0(1), b-fcl = 2,4,... ,n-l. 

To formulate the second theorem we need information on the equilibrium mea- 
sure rf/i^'^^ (x) (see e.g. [SaTop for OP's corresponding to the rescaled weight e~^^"(^) , 
Vn = jjV{cnX + d^), where c^, d^ are the so-called Mhaskar-Rakhmanov-Saff 

(MRS) numbers (sec MhSa, Ra ). For V{x) = K2mX^"' + K2m-ia;^'""^ H as in 

(|1.19|l , we have (DKMVZ2) to any order q as N ^ oo 
(2.3) 

/ 1 (0 "III \ 1/(2™) 1 

\K2Tn m{2m - ly.lj ^ 

and 

(2.4) dN = -^!H!!lzl + y d(,)iV-^/(2™) + 0(iV-(3+i)/(2'")). 

2mn2m ^ 

As iV — > (X), the equilibribum measure is absolutely continuous with respect to 
Lebesgue measure, dfi"^'^^ {x) = V'^'^'' {x) dx, and is supported on the (single) interval 

[-1,1], 

(2.5) #(x) ^ Vjv(x) = ^|1 - a;Y/\hi,i](^) hi,{x) 
(see |DKMV!Z2] ') where hN{x) is a real polynomial of degree 2m — 2 

2m -2 

(2.6) hwix) = J2 ^N,kx^ 

k=0 

and the coefhcients ft.Ar.fc can be expanded to any order in powers of A^^i/t^™) as 
above. In particular, to any order q = 1, 2, • • • , as — s- oo, we find uniformly for x 
in compact sets 

Q 

(2.7) /iw(:e) = ft(a;)+^iV-^/(""U0)(a;) + O(iV-(«+i)/('")) 

where (see lbKMV^2l (2.7)]) 

„2fc . (2m)(2(m-l))...(2(m-fc)) 



(2.8) hix) = J2 f^' 



(2m - l)(2(m - 1) - 1) • • • (2(m - fc) - 1) ' 
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Also by IDKMVZ2I p. 1501, Remark 3] 



(2.9) 



/i(l) = 4m, h{0) = 



4m 



2m- 1 



Note that h{x) depends only on the degree 2to and is independent of the coefRcients 
of V. Note also that h{x) > h{0) > for aU x eR. 



Theorem 2.2. Let V{x) = K2ma;^™ 
J, fc G Z, as N ^ oo 



ie(j)N+j, (t)N+k) 



(2.10) 



CN+AN + ■ {oil), 
0, 



'*2m > as before. Then, for any fixed 



lU-k) + o{l), j - k odd 

j — k even, j ^ k ^ 
j-k = 



where for q — ±1, ±3, 
(2.11) 



we define 



2 . qn cos((7 arcsinx) dx 



In the proof of H2.10(l we compute a double integral over all the different asymp- 
totic regions in K for the OP's {pj} obtained in |DKMV!Z2] (see Subsection 14. 21 et 
seq. below). 

We see from the above, that in a fixed size neighborhood of the diagonal the 
matrices D^o and Coo do not depend asymptotically on lower order terms in V{x). 
Next note that the coefficient in (EHJ satisfies ,DKMVZ2. (2.11)] 



(2.12) bM= [-+oil)]cN, 



aN = dN + o{l) = 0(1), 



N 



Combining our two theorems we can find an explicit expression for Cn, and also 
for all other matrices involved in the (3 = 1 and 4 correction terms in H1.36|l . H1.37|l . 
In particular, as (even) iV — > cxd 
(2.13) 



-D1 9 — 

TV 



1 

2m 



-1 -1 
1 



-1\ 







V -1 -1 ••■ -1/ 

/ I{n) /(n + 2) 
I{n) 

V /(I) /(3) ••• 
/ (o) (^) 



J(2n- 1)\ 




o(l) 



A21 = -D21 



fnH2m 2m-l 
' 22m-l "^JV 



1 



V 



I{n) ) 

■■■ ((«-"l)/2)\ 



o(l) 



The second and the third matrices are Toeplitz. 

By H2.3I) the asymptotic orders of -812, -D21 are exactly opposite. The same is 
true for all the blocks Bij and Dij, i,j = 1,2, and so whenever we have a product 
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of elements from B and A in l|1.36|l . H1.37|l . the product is asymptotically constant. 
Also 

(2m)!! 



(2.14) 



lim 



CjV mK2m 2m-l 



l2m- 



i(2m- 



-1)!! 
^mJ 



1 • 2 ■ 



2(m!) 



M2 



22™-il • 3 • . . . • (2m - 1) 1 • 2 • . . . • m (2m)! 

so that in BA the coefficient K2m cancels out, and as ^ oo, BA converges to a 
constant matrix that depends only on m. 

Recall that in each of the cases /3 = 1 and 4, in order to compute the correction 
terms in H1.36|l . (|1.37|) we need to invert a certain matrix. Now we know that the 
/? = 4 matrix to be inverted, /„ + B12A21, converges as A — > 00. In the P = 1 
case the matrix to be inverted is l2n + BAC, which also converges as A^ ^ 00. It 
is not a priori clear whether the two limiting matrices are related. However our 
structure Theorems 12 . II and 12 . 21 and also Theorem 12 . 51 below, imply the following. 
Let R = Rn denote the n x n matrix with all zero entries apart from ones on the 
anti-diagonal (thus Rij — lifj^n — i + l,l<i<n, and Rij = otherwise). 

Theorem 2.3. yls A^ ^ 00, A even, we have {BA)22 = -R{BA)iiR + o{l). Also, 
quite remarkably, 

/ 

(2.15) BAC = 

V (BA)21+0(1) {BA)22+0{1) 

so that in particular BAC and BA asymptotically have the same two lower blocks. 

It is interesting that the 11 and the 21 entries of BAC are identically zero, not 
just 0(1). 

Now we see that inverting /2„ + BAC boils down to inverting 

/„ - {-R{BA)nR) + 0(1) = Riln + Bi2A2i)R + o(l) 

so that the determinants of the limiting matrices in the (3 = 1 and 4 cases that are 
to be inverted are the same. Thus in both cases we have to check that det(J„ + 
B12A21) = const + 0(1), const 7^ 0. 

Next, we observe that the zero pattern in B12 and A21 in (|2.13|l implies that 

det(/ + B12A21) = (det + 0(1)) • (det T,n-i + o(l)) 

as A^ ^ cxD, A^ even, where 

f 1 1 ... -1\ 

-1 -1 ... 



T 



I'm. 



2(m!) 
(2m) 



1^2 



1 
2m 



V -1 -1 



-1 



V 



(2.16) 



V 



I[n) 
I{n - 2) 

m 



0) VI, 





V 



I{n + 2) 
I{n) 

/(3) 

((n-3)/2) 
1 



I{2n - 1)\ 
I{2n - 3) 

I{n) I 
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and 
(2.17) 



T, 



m— 1 — ^in — l 



2(m 
(2to) 



n2 



1 

2m 



1\ 



1 1 



V 1 1 ... 1/ 

/ /(n) /(n + 2) 
/(n-2) /(n) 

(l) ■■' ((n-3)/2) \ 
(o) "' {(n-5)/2} 

... 1 y 



/(2n- 3)\ 
/(2n- 5) 



V ... 1 

Quite remarkably, we also have 
Theorem 2.4. For any m e N, m > 2, 

detr^ = detTm_i. 

Thus it suffices to show that the single determinant, detTm_i, is nonzero (see 
Theorem 12.61 below "l. The proofs of Theorem 12.31 and 12.41 relv on the following 
identities. Let Dat, denote the N x N matrices formed by the first N rows and 
columns of Doc , Coo , respectively. Recall that by H1.32|l Doc has bandwidth 2n+l. 
Let loc denote the semi-infinite identity matrix. 

Theorem 2.5. (i) The semi-infinite matrices satisfy 

DoC^OC ^QcDoc -^OO 

(the products are well-defined since Doc is banded); 
(ii) the sections for N > n satisfy 



(2.18) 



^nDn — 



i-N-n 








where e^_„ „ is the (N — n) x n matrix formed by the rows 0, ■ ■ ■ , N — n — I and 
columns N, ■ ■ ■ , N -\- n — I of the matrix eoc ■ 

Remark 2.1. Assuming N is even and denoting Wn = £nDn we see from (|2.18f) 
that 



(2.19) = DnWj,\ D-^^ = W^^N. 

Thus the question of the invertibility of both ejsi and -Dat, reduces to checking that 

det(/„ - B12D21) - det(/ + B12A21) ^ 0. 

On the other hand, one can show that the corrections in (|1.36|) . (|1.37|) can be 
rewritten in terms of the lower-right n x n corner of the matrices and D~^ , 
and so (|2.18|) and H2.19|l explain why the same determinant should be checked to 
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be nonzero for both the l3 — 1 and 4 cases. Also using the fact that 

(In-u * 
(In - Bi2D2l)-' 

along with l|2.19|l . the skew symmetricity of eat, and H2.18(l . we can show that e^"'^ 
coincides with Dn everywhere apart from the above mentioned lower-right n x n 
corner. 

Next we state our main technical result that ensures the invertibility of the 
matrices in (|1.37|) in both the (3—1 and 4 cases. 

Theorem 2.6. For any m e N, m > 2, 

detr„_i ^0. 

Remark 2.2. From numerical computations for m = 2, • • • , 14, it seems that det Tm-i 
approaches the value as m grows, but we have not been able to use our methods 
to establish such a result. 

Remark 2.3. The basic idea in proving Theorem 12. 61 is to show that for an appro- 
priate norm || • ||, HTm-i — /m-i|| < 1- Rather than analyzing the quantities I{q) 
in (|2.11() . H2.17|l directly, we consider 

1 2 

(2.20) /((?) = to/(<7)---- sin ^ / coa{qe) y,n{0) d0 , g = 3, 5, • • • 4m-5 
where 

1 1 cos^ ( 



(cf. (|6.3(l et seq. below). It turns out that h{x) is (essentially) a 2F1 hypergeometric 
function (see HCllfl ) and satisfies a first order differential equation (see Ht).8|l '). This 
in turn implies that ym satisfies a Riccati equation 

4 / 2rn + 1 „\ / 1 
(2-22) y'm = —7^ V.n + —— cos 



2 cos 6* 



where we see that the right-hand side conveniently splits into a product with explicit 
factors. Equation (|2.22|l plays a key role in proving that the magnitude of y^ is 
0(^/m), a crucial fact in the proof of Theorem l2.6l In estimating the magnitude of 
the quantities /(g), it is clear from H2.20(l that we are in a double-scaling situation: 
if q is large compared to to, then I{q) is small by virtue of the usual decay of Fourier 
coefficients of real analytic functions. On the other hand if q is small with respect 
to m, the asymptotic behavior of y^ dominates. Both regimes must be analyzed 
separately. These two effects are comparable when q \pm. 



3. Proofs of Theorem II. II and Corollary I1.3I 

We only consider the case (3 = 1. The proofs for (3 = 4 are similar and are left 
to the reader. 

We need the following result which is based on the Plancherel-Rotach type 
asymptotics for OP's in |DKMVZ2] (see Section El for the proof). 
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Theorem 3.1. Let V{x) — H2mx'^™ + ■ ■ ■ as before. Let a = o{cn) as N ^ oo. 
Then, as N oo 

(3.1) UN\\L^a-a,a))^Oic-^'^'), \\eM\ ^M) ^ 0(4^ ■ N-'^^) . 

Also UNh^iw) = 0{c],^^^N^/^). 

3.1. Proof of Theorem [TTTl Set qn = = RnaA'^)- Observe from ((T^ 
and l|4.152|l . I|4. 154(1 that qn satisfies (|4. 13811 with a = 1/2. Thus a* in 14.168(1 is 
1/4. In view of ((T~TT|) . ifT"^ . we have 
(3.2) 

-^SNi(r+ ^,r+ -!^) ^fSNiD)(r + ^,r + -^) 

QN "''-\ ON ' ON/ di,^ "'-^ 'V OJV ' ON/ 



(eS'Ari)(r+ ^,r+ ^) - i sgn(,e - ??) -^Sn i(r + + ^) 

^ On' On/ 2 o '/ ON "'-^V 'on' On/ 



where 
(3.3) 



-Sn.i r H , r H = — Kn r H , r 



(ZJV \ QN QN / QN \ QN QN 

- — (<i>i(r + ^/qNf,0^) ■ (AC(/2„ - BACrY 

QN 

■ (e$i(r + ij/qwf, e<l>2{r + v/qwff 

The convergence of the derivatives and integrals of the ChristofFel-Darboux kernel 
-^Ki^{r + -^,r + ^) with uniform error estimates for £,,ri,r e [—2Lo,2Lq\, for 
a fixed < oo, is established in Subsection 14.31 see 1(4.150(1 . I4.157() . I(4.158() and 
(HTTTli below. 

We now consider the correction terms. By the results of Section|21(see, in partic- 
ular, 1(2.13(1 et seq.), BA and C converge as ^ oo. Moreover J2„ ~ limjv^oo BAC 
is invertible and hence {hn - BAC)~^ = 0(1) as TV ^ oo. From ((^^ and ifT^ . 

(3.4) max |<i>i(r + ^/g^)! = 0(c^^/'), \\e<^2\\L^iw.) = 0{c}^^ N-^'^). 

KM''I<2Lo 

Thus the correction term for the 11 and 22 entries in 1(3.211 is bounded by 

(3.5) const • , , ■ c„'/' • 7Vi-i/(2™) . {c^^N-^'^) ^o(-^ 

as TV — > oo, uniformly for ^,ri,r € [—2Lq, 2Lq\. 

The correction term in the 12 entry of 1(3.2(1 has the form 
1 



1% 



($i(r + a^Nf, 0^) • {AC{l2n ~ BACrY 



■ ($i(r + f]/qNf, ^2{r + v/^N^f 

and is bounded by 

1 

(Arl-l7(2ni))2 

as iV — > oo, uniformly for S^,ri,r <E [~2Lq, 2Lq\ 



(3.6) const . • c-'" • iV-^/(^") • c^/^ = o(l) 
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Finally, to analyze the 21 entry of (|3.2(l (see (|l.;^6(l ). we use the following obser- 
vation. By H1.17|l . {eSNA)ix,y) is skew symmetric. Thus 

(3.7) {eSN,i){x, y) = {eSNA){x, y) - {eSN,i){y, y) = - [ SN,i{t, y) dt. 



After setting x = r + S,/qN, y = r + Tj/qN wc sec that the correction term has the 
form 

i-r+ri/qN , T \ 



(3.8) Jr+ilqf, 



(e$i(r + 7?/g7v)^, £$2(7- + v/flN?^^^ 



The integration is over a subset of the region (see (|4.6f) below) . Hence as in (|3.5|l , 
but now using H4.18|l below, H3.8|) is bounded by 

(3.9) const • 4^N-^ ■ N^-^/^^"'^ ■ c]I^N"'^^ = 0{N-^'^) 

as iV ^ 00. We compute finally the contribution of the ChristofFel-Darboux to 
(|T7|) . We have 

(3.10) - r KN{t,y)dt ^ {eKN){x,y) - (eKN){y,y) 

where x ^ r + £,/qN, y ^ r + tj/qn. By H4. 159(1 . 14.171(1 below with a = 1/2, a* = 
1/4, 1(3.10(1 equals as ^ 00 

(^1^ \-^{t)dt~£ \^(t)dt^ +0{N-^/^)^ 1^ \^{t)dt + 0{N-^/^) 

where the estimates are uniform for r,^,ri£ [— 2Lo, 2Lo]. 

Assembling the above estimates we see that we have proved Theorem II .11 in the 
case /3 = 1 with error term 



(3.11) £n.i 



as N 00, uniformly for £^,ri,r E [—2Lo, 2Lo]. 



3.2. Proof of Corollarv ll.3[ In (TW2j the authors prove that if xb denotes the 
characteristic function of a bounded Borel set B then for finite A^ in the case (3—1 



Prob{ no eigenvalues in B } 



\ 



det (1 - Kn,iXb 




where again e denotes the operator with kernel i sgn(a; — y) . Now the product of the 
operators above is of the type identity plus trace class, but the individual operators 
are not. This is because e is Hilbert-Schmidt but not trace class. However, as 
indicated in Section ^ we can use regularized 2-determinants. One easily sees. 
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moreover, that det2 



det {1-Kn,iXb) 



1 






^ 1. Thus 










( ' 










det2(l - Kn,iXb) det2 _ 
For /3 = 4 such regularization issues do not arise. For P — \ (cf. (11.301) 1 

(3.12) 



Prob < no eigenvalues in , , 2 ' ' ' \ 2 



For P — 4, the same formula is true provided we replace Aat^i Aat^, det*-^-* 
det^''^ Kn,i KNf2,4 (cf. III23). In order to prove Corollary 1 1.31 for /? = 1 it is 
clearly sufficient to prove that the entries of ^iC^"'^^ (r + r + = 
apart from the term i sgn(^ — 77), converge to the corresponding entries in K^^\^, 77) 
(cf. lEiSl, EiSl), C^) in trace norm for L'^{-0,9). The case /3 = 4 is similar, 
mutatis mutandis, and is left to the reader. 

We consider r, ^, 77 in the compact set [— 2io, 2io] with Lq = 9. We have 

1 C ^ -i ^ ^ 
— S'at 1 (?' H , r -\ ) — — (r H , r H ) + correction term. 

Qn Qn Qn Qn Qn qN 

By (|3.3|l . the correction term corresponds to a finite rank operator, and hence its 
trace norm is bounded by the L^{—9, 9) norms of the vectors (<i>i(r + ^/qn)'^ , 0"^), 
(e<i>i(r + Tj / q^)'^ , e$2(?' + v/in)'^)'^ which in turn may be estimated by the bounds 
in l|3.4|l . We conclude that the trace norm of the correction term decays as N~^^'^ 
uniformly for \r\ < 29. 

The fact that -^Kn^t + + ^) ^ ^00 (C ^ v) in trace norm was proved in 
the case V{x) = x^™ in [D]. We now give a proof for general V{x) — H2mX^™' + • • • , 
which in addition also applies to the other entries of ^fC^"'^''. Note first that 
it is sufficient to prove that x-^K^f^r + —,r + converges to xKooi: — ■■)x 

in trace norm for any C^{— 29,29) function x with xiO = 1 for ^ G (—9,9). Let 
P = iD = i-^ denote the self-adjoint operator of differentiation with periodic 
boundary conditions on [—29, 29]. Then D + I = \P + / is an invertible operator 
with eigenvalues 7^ = |^ + 1, fc e Z. In particular is Hilbert-Schmidt on 
L'^{~29,29). Clearly x^KnX maps L'^{~29,29) into the domain of P, DomP = 
{/ e L2(-20, 261) : /' G L2(_26I, 29), f{-29) = f{29)} and so 
(3.13) 

X—KnX - JT^iD + I)x—KnX 
QN D + I qN 

x—KNX+-Fr--px'—KNX + -Fr--px(-K7—^N]x 



D + I^qN ^ D + I^qN D + I^\dUN 
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But by (|4.15()|l below, -^Kn + ^ ' + ^) converges to Koo {£, ~ i]) uniformly for 
k|j ICIi 1^1 ^ 20 and hence x^ifjvx x^ooX iu Hilbert-Schmidt norm. But then 
D+j X^-^wX ~* "DqTf X^ooX in trace norm. The same is clearly true for the second 
term in and also for the third term (take j = 1, fc = in (I4.15U|) below). This 

shows that -^K^ — > K^o in trace norm and completes the proof that the 11 entry of 

^i^AT (r+ ^) ^ Koci^, — '']) converges in trace norm to {K'^^^)ii, as desired. 

Clearly the same is true for the 22 entry. Similar considerations using H4.157|) in 
place of H4. 150(1 prove the corresponding result for the 12 entry. Finally, for the 21 
entry we observe that Yn{£^, 77) in (|4.171l) is a sum of two terms Aj^{£_, 77) + Bj^{ri). 
The first term can be controlled in trace norm using the operator P — iD as above, 
together with the estimate (|4.171f) fa). On the other hand B^iv) corresponds to a 
rank one operator in L^(—9, 0) and hence its trace norm goes to zero as — > c» by 
(|4.171|l fb). Finally, the correction term in the 21 entry can be controlled as above. 
This completes the proof of Corollary ^31 

4. asymptotics of integrated op's, proofs of theorem 12.11 and 

Theorems [O IXTI 

Notational remark: Throughout this Section, the notation N, P^oo means 
that N,P ^ 00, N > P, N - P is fixed (and finite). 

4.1. Proof of Theorem 12.11 Let V{x) — K2,nX^"^ + • ■ ■ , {Pj}j>o, {4'j}j>o be as 
before. The following observation is due to A. Borodin. 

Lemma 4.1. For any j,k > 

(4.1) (i?0j>fe) = I sgn(j - k) {V'^„<l>k). 
Proof. For j < k we have 

{Dcj^j,M^ I (p,(x)e-^(-)/^)'p.e-^(-)/2dx 

= J p; (x)pfc(x)e-^(") dx - i(l/'0„ 0fe) 

and the last integral vanishes since pk is orthogonal to all polynomials of degree 
< k, and we have assumed that j < k. Observe that (F'(/)/c, 0fe) = as it should 
be. The case j > k follows by skew symmetry. □ 

Remark 4.1. Clearly Lemma [4. II is true for any differentiable V{x) such that 
J (1 + \V'{x)\) e-^(^' dx<oo, g = 0, 1, 2, • • • . 

As V'{x) is a polynomial of degree n — 2m— 1, it follows immediately from ((4.1|l 
that {{D(j)j,(j)k))j,k>o is a banded batrix, {D(j>j,(f)k) = for \j — fc| > n. 

For the proof of Theorem 12.11 we need to evaluate the leading asymptotics of 
iy' 4)N+j,4'N+k), j,k fixed, as A^ ^ 00. This clearly reduces to computing the 
large A^ asymptotics of {x''4>N+j,4'N+k), j, k fixed, < q < n. In order to compute 
such asymptotics we will use the following relations for the recurrence coefficients 
in (E3J 

(4.2) bN/bN-i^l, and qn — o{bN), as A^ ^ 00. 
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These relations certainly hold for our potential V{x) = K2mX^"^ + • • • , see H2.12II . 

(lOl . 1231). 

Lemma 4.2. For any fixed q G N, as N ^ oo 



x'^(t)N{x) =b% 



(4.3) 



+ Eo(l) • (t>N-q+{2l-l)ix) 



1 = 1 



where the notation o(l) indicates terms which are independent of x and uniform 
for < I < q. (The o(l) terms above can be replaced with 0(iV"^/''^'"^).j 

Proof. By (|4.2|l the three term recurrence relation (|2.1I) takes the form as iV ^ oo, 



(4.4) 



X(f>N{x) = Bn [(f)N+lix) + (1 + oil)) (l)N^i{x) 

+ o{1)(I)n(x)] 



which coincides with H4.3|l for q — I. A simple induction on q using (|4.4|l and 
Pascal's triangle, now gives the result. □ 

We are now in a position to prove Theorem 12.11 Write V{x) — ^ 
(4.5) 



2m 

Then by Lemma [4. II Lemma [4.21 and the orthogonality of the (pj^s, 



^ 2m 

{D(l)N+j,(l)N+k) =2 sgn(i - k)^qVq{x'i^^(l)N+j,(l)N+-k) 

9=1 

^ 2m 

= :jsgn(j - k)^qvq0}^ 

^ 9=1 



/=o 



(l + o(l)),52;,fe- 



-j+q-l 



9-1 

1 = 1 

as — > oo . Observe that for | j — A: | < n and q — 2m = n+1, k — j + q — l — k — j + n 
lies between and 2n. Hence ii \ j — k\ < n and j — fc is odd, we see that the square 
bracket in (|4.5() for q — 2m gives rise to a leading contribution {(n+k-j)/2} + 
However if \j — k\ < n and j — fc is even, the contribution is o(l). But by H2.12II . 
(|2.3|) . — ^ oo, and hence the terms corresponding to g = 1, 2, • • • , 2m — 1 in (|4.5|) 
contribute to lower order. This completes the proof of Theorem l2.1l 

4.2. Integrating the Plancherel Rotach type asymptotics for OP's and 
the proof of Theorems IsTTI and [2T2I 

4.2.1. Auxiliary estimates on the integrated Plancherel-Rotach asymptotics. For the 
convenience of the reader, we recall relevant results from D KMVZ2| . Assume that 

V{x) = K2mx'^"' H , m g N, K2m > 0. Let pn be the Ath OP on K with respect 

to the weight e-^(^) and set (/)jv(a;) = pN{x)e~^^'''>^'^ as before. In |DKMV2!^I 
Theorem 2.2] an asymptotic expansion, as A — > oo, is derived for (jj^ in the whole 
complex plane, and the leading term in the expansion in each region (see below) 
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is computed explicitly. Here we are only interested in real values of the argument. 
Fix So>0 and sufficiently small (see jDKMVZ2j 'l. For < S < So define 

/5 = (-c^,-I-(5), h=i-l-5,-l + S) 

(4.6) h = i-l + S,l-5), I2 = {l-S,l + S) 

/i = (I + S, +oo) 

and let hiq{x), h{x) be as in (ESI, (I22I) and |(2!Hll, respectively. In |DKMVZ2j . 

each of the regions I2 and is split further into two regions, I2 = Ci^s U C2,5, 
I4 — Di s U D2.S, the estimates in C2,s and 1)2.(5 being finer than those in Ci^s 
and Di^s, respectively. For our purposes it is sufficient to conflate the regions and 
use the estimates in Ci^s and Di^s for all points in I2 and I^, respectively. The 
apparent singularities in the formulae below for z — ±1 are of course removable 
(see e.g. (|4.1()|l and H4.12(l below). In what follows, notation of the type 0{1/N) 
as TV — !■ cx), means that the estimate holds for x uniformly in the respective region 
(|4.6|) under consideration. Also the estimates are uniform for 6 in compact subsets 
of (0, Sq]: in this connection see the important Remark |4 . 31 following H4.136() below. 
Finally, the constants ci,C2, ■ • • below may depend on 6, but they are independent 

of X and N, P. 

We are ready to state jDKMVZ2l Theorem 2.2]. First, we consider the "expo- 
nential" regions /i.s. Uniformly for all |a;| > 1 + (5 

4'n{cnx + (In) =- 



V^TTCN 

X- 1 



1 



1/4 



1 



1/4 



1 + 



a; > 1 + ^ 



and 



(j)N{cNX + (In) 



•s/47rcjv 
V x-i 



1 



1 



1/4 



1 



1 



X < -1 - S. 



X e 



N > Ni{V) 



In particular, as 

(4.7) hN{x) > > 0, 

r |DKMVZ2[ Prop. 5.3]), we have for \x\ > 1 + 5 



(4.8) 



\4>NicNX + dN)\ < 



v/47rcjv 

where we have used the fact that jy^ — > C2\y\ for > 1 



Next, consider the 
1 + S 
(4.9) 

4'n{cnX + ^Af) = 



where 
(4.10) 



■'Airy" region on the right, 12- Uniformly for al\l — 5 < x < 



X 



1 



X - 1 
X - 



1/4 



\fN{x)\^'^kiUN{x))(l + 



X - 



1/4 



l/iv(x)|l/4 

fN{x)=aNN^'^ {x-l)fN{,x) 



M'{fN{x)){l + 
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and in turn (see (the proof of) IDKMVZ21 Proposition 7.3]) 

(1) fnix) is real analytic on (1 — 5, 1 + 5), and to any order g = 0, 1, 2, • • • 

fN{x) = ^iV-^/(2™) /(^.)(x) + 0(Ar-(«+i)/(2™)) 
i=o 

uniformly for x in the interval. Moreover, the functions /(j) (x) are also real 
analytic oii\~5<x<\ + 5 

(2) to any order q = 1, 2, • • • 

aN = (/i^(l)/2)'/' = 2m2/3 + ^ iV-J/(2™) ^^^.^ + o(7V-(«+i)/(2'")) 

i=i 

(3) /]v(a;) = -a^N'^/'^UNix), where C/iv(a;) = fN{x) + (x - l)f'j^{x) also has 
an expansion uniform in x to any order (7 = 0,1,2,-- - as above 

q 

Un{x) = ^ Ar-J'/(2™) [/^^.^(a;) + o(A^-(«+i)/(2r«))_ 
i=o 

The terms U(^j) (x) are real analytic onl — (5<a;<l + (5 

(4) maxfc=o,i,2maxi_5<:r<i+5 |d''7Af(2;)/c^a^'' I < M < cx3 for N > N2{V) 

(5) /jv(l) = 1 = Un{1) and mmi^s<x<i+s fN{x) > i for TV > Ar2(V^). Also 
/(o)(l) = l = C/(o)(l). 

Similarly, in the left "Airy" region I4 we have uniformly for all— 1 — (5<x<— 1 + 5 

(4.11) 

cPN{cNX + dN) =^-^(^\^\'^'\fN{x)\'/'Ai{-Mx)) (1 + 



X + 1 



X- 1 



' AiVM-^))(l + 0(l))) 



l/iv(x)|l/4 



where 

(4.12) /Ar(x) = (5a,7V2/3 (x + 1)/Ar(x) 

and in turn 

(1) fN{x) is real analytic on (—1 — 5, — 1 + 5), and to any order q = 0, 1, 2, • • • 

/jv(x) = ^iV-J/(2™) /(^.^(:r) + 0(Ar-(9+i)/(2™)) 
i=o 

uniformly for x in the interval. Moreover, the functions /(^-j (x) are also real 
analytic on— 1 — (5<x<— 1 + 5 

(2) to any order q = 1, 2, • • • 

aN = {hl{~l)/2f'^ = 2m2/3 + ^ iV-J/(2™) + ©(Ty-^^+D/^^™)) 
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(3) fj^{x) = -aNN^^^UNix), where UNix) = /^(x) + (x + l)/^(x) also has 
an expansion uniform in x to any order q = 0, 1, 2, • • • as above 

q 

i=o 

The terms U^j) (x) are real analytic onl — (5<x<l + (5 

(4) maxfe=o,i,2max_i_5<.^<_i+5 |d'=/Ar(x)/dx'= | < M < oo for N > N2{V) 

(5) = 1 = UNi-1) and min_i_,-<,<_i+5 /Ar(x) > i for iV > N2{V). 

Also/(o)(-l) = l = C/(o)(-l). 
Finally, in the middle region I3 we have uniformly for all— l + (5<x<l — 5 
(4.13) 

(I)n{cnx + dN) = \ -Tj JTY74 

V ncN |1 — x^|-^/* 

X (cos(y^ |l-2/2|i/2/,^(y)rfy+ ^arcsinx) + 

+ sin(y ^ |l-2/'r/'/iJv(2/)dy- ^arcsinx) 

In the proof of Thcorem l2.2l fand Thcorcm l3.1|) we need the following two results 
that hold as — > 00, 
(4.14) 

4>N{y)dy = CN I 4>N{cNy + dN)dy 



and uniformly for x G K. 



(4.15) 



CN / (f>N{cNy + dN)dy 



= 0(cf A-i/2). 



or, equivalently, for any interval K dM. 

CN / 4>N{cNy + dN)dy 

JK 

Notational remark: In (|4.14() above the reader may wonder why we write 0{N~'^/'^)+ 
0(iV-i/(2m)) rather than just 0(A^-i/(2™)). The reason is that the 0(iV-i/(2'")) 
term arises simply from the evaluation of the constants a at, oin arising in the asymp- 
totics using ^^{2), lin^ f2). The term Oi^N''^/'^), however, constitutes the de- 
tailed estimate of errors involved in evaluating the integral over all the different 
asymptotic regions (|4.6|l . In order to separate out these two very different sources 
of error we have adopted the convention of writing error estimates in the form 
(9(jV-i/2^ _l_ 0(7V-i/(2»™)) as above. We use this convention in all the estimates 
that follow in this Section. In the proof of Theorem 12. 21 error estimates of the type 
0(iV^^/(2m)) also arise from the approximation of /lAr(x) in (|2.7() (see e.g. H4.129() 
below) . 

Note that (|4.14|l . H4.15|l are direct consequences of the following more detailed 
statement, which is also used in the proof of Theorem 12.21 
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Proposition 4.3. The following holds uniformly for x in the respective regions as 
N oo. In Ii: 



(4.16) 



in h •■ 
(4.17) 



CN I 4iN{cNy + dN)dy 

1+5 



<CN \(t)N{cNy + dN)\dy 
Ji+s 

= 0{{cN/N)^/^e-^^'), x>l + S; 



i+s 



CN / <t>NicNy + d^) dy = c]{^N-'/\2m)-'^^{l + 0{N-^'^) + 0(^-i/(2"))) 
Ji-s 



CN I 4>N{cNy + dN)dy 

1-5 



< o 



Cn 

Ari/2 



in h 



(4.18) 



CN I 4>NicNy + dN) dy 
i+S 



l-6<x<l + d; 



,1/2 X 



<0(^j, -l + S<x<l~S; 



in I^: 
(4.19) 

CN f ' ^NicNy + dN) dy = 4/'7V-i/2(2m)-i/2(-l)^(l + 0{N-'^') + 0{N-^''^^'' 

J -1-5 



-1+5 



CN I (f)NicNy + dN)dy 

1-5 



< o 



-l-5<x<-l + S: 



in A ; 



(4.20) 



Cat / (t)N{cNy + dN)dy 



< CN \4>NicNy + dN)\dy 

J — oo 

= 0((cAr/iV)i/2e-^=^), x<-l-5. 



Proof. We refer to (|4.8|) and note that 



1+5 



which proves (|4.1t)|) . The proof of H4.20|l is similar. 

We now prove H4.18|l . Uniformly for — 1 + i5 < a; < 1 — (5, by (|4.13|l . 



„l/2 

-N 



■ arcsm y 



CN / (t>NicNy + dN)dy 

J -1-5 

[\^-i'\"'hN{t)dt+\ 
+ sin |1 - i'\''''hN{t)dt - i arcsiny) o(^) 

/ N n 1 

\2 J ~ t'\''^hN{t)dt + - arcsiny 



dy 



.1/2 

-N 



1+5 



|l-y2|l/4 

dy 



|l_y2|l/4 



+ 0(cfiV-). 
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Integrating by parts in the last integral we find 

sin 



-N 



arcsm y 



-1+5 



dy 



|l-y2|l/4 

2 arcsiny] 



N 



|l_y2|3/4/j^(y)+ l|l_y2|-l/4 



(4.21) 



y=-l+5 



.1/2 

-N 



sm 



1+5 



\2 J \^^^''\^'^^N(t)dt+-'Avcsmy 



dj/Vf |l-y2|3/4/j^(y)+l|l_y2|-l/4 



dy 



uniformly for x £ I^, which proves H4.18|l . (We have used 14.7|l and also the uniform 
boundedness of h'j^{x) on /a, as TV ^ cx).) 

Now we prove (|4.19(l . Uniformly for -1 - (5 < .t < — 1 + (5, by (|4.11|l . 



CN 



(j)N{cNV + dN)dy 



(-l)V/' 



1-5 



y-1 



y + 1 



1/4, ~ 



\j^{y)\^/^M{-f^{y))[l + 0[- 



Now 





y+l 


1/4 




y- 


- 1 








y- 


- 1 






y- 


h 1 



1 



\fN{y)V'' 

1/4, ~ 



■Ai'(-/^(y))(l + o(^)) 



dy. 



\}^{y)\^/^M{-h{y))0[^)dy 





y - 


1 


^ J-1-5 


y + 


1 



1/4 



dy = Oi4^N-') 



where (cf. HESt]) 



CiEEsup|2/r/4|Ai(y)| <oo. 



Next, in view of H4.12|l and the properties of /jv, 



,1/2 



1-5 



y + 1 



i/4AiVMy)) 



O 



/ 1 , l+i 



< const ■ CjI'^N 



5 
1+5 



-1-i 
-H 



y+1 



y-1 

y+l 



l/iv(2/)|i/4 
l/4 1 + |a^jV2/3(y + l)/^(y)|l/4 
|5jvW2/3(y+i)/^(y)|l/4 



dy 



^ ^ (7V-V6 + |y+i|l/4)^y^O(,V2^-l)^ 



where (cf. HESi]) 

(4.22) C2^sup(l + |y|i/4)-i|Ai'(y)| < 
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and we have used (see jDKMVZ2l (8.72)] and lfTJ^ (2)) 

(4.23) UN = (/i?^(-l)/2)'/^ = 2m2/3 + 0(iV-i/(2™)), 



We consider next 



(4.24) 



1/2 



dy. 



l-i-s y ~ 1 
Make the change of variables (cf. I|4.12|l . (|4.23(l 'l 

(4.25) uiy) = -fNiy) = -5ivA^'/'(2/ + l)/iv(2/) 

which iniphes (cf. iflTT^ fS')') 



(4.26) 

Note also that 
(4.27) 

Now set 
(4.28) 



du 
dy 



UN{y) = fN{y) + {y + i)fN{y). 

— 7/ 

y{u) + 1 



BNiy) 



1 



\y-i\'^VN{y)\'^'UN{y) 

where the denominator is bounded away from by (|4.12|) for 5 > small enough. 
Then (|4.24|) becomes 
(4.29) 



.1/2 M^) 



-AT 



Ai'(ii) BN{y{u)) du 



1/2 



k\{u)BN{y{u)) 



u{x) 



, B'{y{u)) , 
Ai{u) . " du 



u(-i-5) o^nN^/^ Jui-i-5) UN(,y{u)) 



since the boundary term and the last integrand are uniformly bounded for x I4 
and the length of the interval of integration is 0{N'^^^). 
Finally we consider the integral of Ai. Set 



(4.30) 



\y-i\'^V~N(y)\'/' 

Un{v) 



and for future reference note that (see H4.12ll f5)) 
(4.31) Ln{-1)^2^/\ 
Making the same change of variables (|4.25(l , H4.26(l , l|4.27|l we find 



(4.32) 





y + 


1 


1/4 


1-1-6 


y - 


1 





\fM{y)\"'M[-fr,{y))dy 
Ai(u) Lj\j{y{u)) du. 
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In what follows we use the standard estimates (see jAbStp 

I Ai(t)| < Ce"(2/3)t^/^^ I < cil + tl/4)g-(2/3)t^/^^ ^ > 0, 



(4.33) 

and 
(4.34) 



|Ai(i)| <c(l + |^|)-l/^ |Ai'(t)| <c(l + |^|)l/^ i<o 



Ai(s)ds < c(l + |^|)-^/^ t<Q. 



Note that by l|4.25|l . u{y) > for y < — 1 and u{y) < for ?/ > —1. Consider first 
for X < ~\: 



(4.35) 



.1/2 



(-1) 



< 



N 



m(-1-<5) 



[x) 



Ai(u) Lisi{y{u)) du 



„i/2 

-TV 



4/^7VV2 



const 



\A{{u)\du^O{cjl^N-^'^) 



Next we consider a: ~ —1. Note that m(— 1) = and \m\M^oou{—\ — 5) ~ +cx). 
Also for any fixed u e [0,u(— 1 — 5)], hniAr^oo -^jv (?/(«)) = Li^{—\) ~ 2^/'' (see 
llOTIl . E3ni), 113U, and also llO^ V We have 



u{-l-S) 



Ai(u) LN[y[u)) du = Li^{-1) 

ru(-l-<5) 



n(-l-(5) 



Ai(u) dit 

Ai(w) {LN{y{u)) - Ln{-1)) du = Qa,i + Qo,: 



Now by K77\i 



\Qo,2\ < 



|Ai(u)|(_^ 



< c 



u(-l-(5) 



|Ai(«)| 



max \L'^it)\]\yiu) + l\du 

-o<r< — 1 



\aNNV3f{y{u))\ 
<cN-^/^ I \uAi{u)\du = 0{N-^/^). 



Also Qo4 = 21/'*(/q°°-£[1^__5j) Ai(u)(iw. Taking ll05ll . (lO^ . ^J^{2) into 

account, we conclude that 

(4.36) 

l/2^_-|NAr „„(-l_5) 



Ai(u) Ljv(y(u)) 
=1/^ 21/4 



Now consider —l<x<—l + S. We have u(a;) < 0. Then 14.32|l becomes 



(4.37) 



Ai{u) Lj\i{y{u)) du 
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where the first integral was evaluated in (|4.36|l . and for the second we write 
(4.38) 



u{x) 



PJi{u) Lpf{y{u)) du 



— oo 

, nu 



Ai{t)dt) LN{y{u)) 



u{x) 



■ du 



u{x) ' ^ — oo 




o 



du 



(1 + 1^1)3/4 iV2/3 



by H4.34|l and the uniform boundness of L'pj. From this we draw two conclusions. 
First, uniformly for —l<x<—l + 6 



(4.39) 



Ai{u) LN{y{u))du = 0(1), N oo, 



and, second, again using (j4.34|l . for any fixed —l<x<—l + d, and in particular 

for X ^ —1 + S, 



(4.40) 



A: 



u{x) 



(u) LN{y{u)) du = 2i/4( / Ai(t)dt) + 0(iV-i/2) 



(4.41) 



i{x) = -aNN'^^^ix + 1)/Ar(x) < -C{x) N'^/^ -OO. 



Hence we conclude that the uniform estimate in (|4.19|l holds. Also recalling H4.36|l . 
(O^ we find 



'^N \ ^) r)l/4 I 



r Ai{t)dt+ f Ai{t) dt + 0{N-^^^)'^ 



1/2 

Ari/2 



since Ai(t) = 1 (see |AbSt| ). which proves the asymptotic formula in H4.19|l . 
The proof of (|4.17|) is similar. □ 

4.2.2. Proof of Theorem \H.l\ First we prove the second statement. By the definition 
of e 



(4.42) 



1 r r 

e(j)N{x) ^ - / (j)N{y)dy 



+ C30 



(t>N{y) dy 



2 r+oQ p+oc 

/ (j)N{v)dy- I 4)N{y)dy. 



Making the change of variables y c^y + ^at, we see that we have to estimate 
(4.43) 



-CN I (f>N{cNy + dN)dy 
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and 

(4.44) CN / 4>N{cNy + dN)dy 

J {x~dN)/cN 

The second statement now follows from ()4.14(l . H4.15|l above. 

Now set y — cnx + i-e. x = {y — dN)/cN. Note that if ?/ S [—a, a] where 
a = o{cn), then x — o(l), and hence 4'N[y) = 4'N{cNX-\-dN) is given asymptotically 
by (|4.13|) . Thus \\4>N{y)\\L°^([-a.a\) — 0{cj^^'^) which proves the first statement. 

Finally to prove the third statement (which is not used in this paper, but is of 
independent interest) we again introduce for y the rescaled variable x,y = CNX+dN, 
and note that by (|4.8|l . 

sup \(I)n{cnx + dN)\ = 0{c^''''^ e'^"^), 

\x\>l+5 

and by PTl|l . 

— 1/2 

sup \4>N{cNX + dN)\^0{c^' ). 
\x\<l~S 

Thus ||'/>jv||l°=(r) is determined by the Airy regions. Let us consider —1 — 5 < x < 
— 1 + 5 (the neighborhood of x = 1 is treated in the same way). We refer to (|4.11() 
and note that by (|4.25(l the Ai term is of order 
(4.45) 

\c~''^\x + l\-''^\aNN^'\x + 1)In{x)\^'' ki{u{x)){l + 0{N-')\ = 0(c^'/'iVi/6) 

since Ai is bounded (and is not of a smaller order since e.g. Ai(u(— 1)) = Ai(0) ^ 0). 
Also by 14.22|l for the Ai' term 

^V2|^_^i|i/4AiVM£0)(i + o(^- 

^ |/^(:e)|i/4 

c.c-/^|. + i|i/4 (i + l/;v(x)r) ^ ^ ^^^-^1/2) 



(4.46) 



~ ^ l/iv(x)|V4 



by (05l) . 

Remark 4.2. Note that in [Sz] the estimate ||0Af||L°°(M) — 0{N^^/^^) is proved for 
the case of Hermite polynomials, V{x) — x^. Of course in view of (|2.3(l for m — 1 
this is consistent with our estimate 0(c^^^^iV^/^). Also from our estimate and (|2.3() 
it follows that, for all m > 2, ||0jv|1l°=(r) blows up as Afi/6-i/(4ni)^ as TV ^ oo. 

4.2.3. Proof of Theorem \2.S\ From l|2.10|l we anticipate the leading order terms to 
be of order ^ - ^j0^- Using (jO^ we find 

{(I^N , e(j)p) ^^(^ J <j)N{x)dx^(^J <j)p{y)dy^ 



(4.47) 



— OO 

+ 00 /'+00 

(j)Nix)dx / (t)p{y)dy. 



OO 



By (|4.14|l . as N, P^oo, the first term in H4.47|l equals 

(4.48) 2 {NPy/^ 2m ' ^ J + ^l^v )) 

X (1 + (-1)^ + 0(p-i/2) _^ Cl(p-i/(2™)))_ 
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The estimation of the second term in H4.47|l involves an extended region-by-region 
calculation, and concludes eventually at equation H4.136|l below. 

The following result will be used repeatedly throughout this Subsection. 

Lemma 4.4. For V{x) — H2rnx'^^ + ■ ■ • > i^2m > as before, as N, P^oo 
(4.49) 

and 



cp 2m P 



(4.50) 



oJjv — dp 

Cp 



0{P 



-l-l/m 



Proof. Note first that for any fixed a > 0, as N, P^oo 

1 1 _ 1 11 

(4.51) TV" ~ ~ {P+{N - P))" ~ 



N-P ^ 



1 

P2 



oip-'-") 



This together with (|2.3|) for q ~ 2m implies 



cn 



CP =C(_i) [7Vl/(2™) _ p(l/2m)] _^ J2 [C(fe) • (TV-'^/^^") - p-fe/(2rn))] ^ 0(p-l-l 



/(2 



fc=l 



C(_i) [iVl/(2™) _ p(l/2m)] _,_ J2 [C(fe) • 0(Ar-l-'=/(2™))] ^ Q(p-l-l/(2™))^ 



fe=l 



Therefore by (|4.51|l 

CN - CP _C(_i) [Arl/(2m) _ p(l/2m)] + ^(p- 1- l/(2m) ) 



Cp 



C(_i)Pl/(2rn)(l + 0(P-l/(2m)) 
[((iV/P)l/(^™) - 1) + 0(P-1-1/™)] (1 + 0(p-l/(2m))) 



1 AT-P 



2m P 
1 iV-P 



0(p-2) +0(p 



-l-l/r?x\ 



(1 + 0(P 



-l/(2m) 



)) 



+ o(p^-^^^/^^™^). 



2m P 

The second statement follows similarly from (|4.51|) and (|2.4f) with q = 2m. □ 

We begin the estimation of the second term in (|4.47l) by making the change of 
variables x cnx + (In, y cpy + (In in 14.4711 to find 
(4.52) 



+00 



(t)N{x)dx I (t>p{y)dy = CNCp j (j)N{cNX + dN) dx (j)p{cpy + dp) dy 

) JX(x) 



+ 00 



CNCp 



^ /> /'+00 

y2 / (I)n{cnx + dN)dx / (j)p{cpy + dp) dy 
.=1 •'x(^) 



where 
(4.53) 

Note that by Lemma lOI 



Cn dN — dp 
X{x) = X ! . 

Cp Cp 



Cn _ ^ 

Cp 



N- P 1 
2m P 



^0(p-l-l/(2m)) > 1, 



dN — dp 
Cp 



= o{p- 



-l-l/r. 
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and so for any fixed xo > (e.g. xq = 1/2) and N, P large enough 

(4.54) ^{x) > ^ for a; > xq, ^i^) < x for a; < —xq. 

Now gjH), g^OJ, imply 
(4.55) 

i + )4>N{cNX + dN)dx / (j)p{cpy + dp) dy 

^J/i Jl.y Jx(x) 



CNCp 



< CNCp 



] ■ sup 


L 


/ xeK 





r'+oo 

cj)p{cpy + dp)dy 



{ i l'^ )\4>n{cnx + dN)\dx) ■ sup 
= 0((cwcp)i/2(iVP)-i/2 
Hence in H4.52|l we are left with three types of integrals 

CAfCp / + / + / )(?!'Ar(cAra; + rfw) do; / (l)p{cpy + dp) dy 
(4.56) "'/s Jh/ Jx{x) 

= J4 + J3 + J2. 

By H4.54|l . uniformly for x G I4, we have X{x) < x for N, P sufficiently large. So 
(4.57) 

J4 = CNCp / (/)Ar(cAra; + dAf) da; 



/ +[+[+[ )(l)p{cpy + dp)dy 

Jx{x) J la JI2 Jh J 



Jh 

= J44 + J43 + J42 + J41- 

Now J43, J42, J41 are products of (|4.19(l with H4.18;i . (|4.17(l . H4.16|l respectively, and 
hence 

(4.58) J42 = ^'2^'' (1 + + 0(7V-i/(^"))) 
and 

(4.59) J43 + J4i = ^_-^0(P / +e )^___0(P /). 
Next define x^^p G /4 such that 

C/V C?A/^ — (ip 

(4.60) X{xNp) = XNP—^— - = -\-6. 

cp cp 

By Lemma [4.41 we see for iV, P^oo 

(4.61) x^,p = -(1 + S)(l- + 0(P-1-V(2™)) 

and so xn p indeed lies in I4 = (—1 — (5, — 1 + S). With this notation 
(4.62) 



J A. 



xn.p /--l+iS , /--l+iS 

cn4>n{cnx + dN)dxi / cp(j)p{cpy + dp)dy\ dx 



— J44 + J44. 

Observe that in J44, —1 — 5 < X{x) < —1 + S and so all the points y in the inner 
integral lie in I4. 

We need the bounds 

1 1 1/2 

(4.63) sup \4>n{cnx + dN)\ < C Cj^ ' 

min{\x+l\,\x-l\)>S 
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and 
(4.64) 



C 



which were almost proved in H4.45|l . H4.46|l . Indeed, the two inequahties preceding 
imply (|4.t)3|l . and as 



I Ai(/^(X))| < C • (1 + \fN{x)\)-^'\ I Ai'(/Ar(x))| < C • (1 + \fN{x)\Y'\ 



we have for |a; — 1| < 5, by H4.10|l . 
(4.65) 



< 



C|/Ar(a:)|i/4 



< 



\f^{x)\^/^\Ai{f^{x))\ 
l|i/4 

\x - l|i/4| Ai'(/jv(x))| ^ C(l + |x-l|iV^/3)V4 

< r^T^ < const. 



C(Ar2/3)i/4 



Ix - l|i/4(i + |/jv(a;)|)i/4 - (1 + - l|iV2/3)i/4 



Ari/6 



Inserting H4.65|l and the analogous inequalities for |a; + 1| < S into 1)4.91) . we obtain 

gSl- 

Now by Unm . 

^44 < Cat / 777; ,^ dx 



-1-5 C 



1/2 pl/2 



const ■ 



But by l|OT|l 

Thus 

(4.66) 

Next 
(4.67) 



, 1 , X dp - (In cp - cn ^/n-u 



cp 



cp 



W < const . A.Vap-1 ^ ^(^-1/3). 



4'4 = (-l)^C 



AT 1/2 



-1+5 



x-1 



x + 1 



1/4, ~ 



|/^(x)r/4Ai(-/Ar(x))(l + 0(l 



a; + 1 



x-1 

-l+S 

cp<i)p{cpy + dp)dy 

X(x) 

We consider first the two 0{N^^) terms. As 

(4.68) \\fN{x)\"^Ai{~h{x))\<C' 



1/4 1 ~ / / 1 

^ -Ai'(-/w(x)) (l + 0( — 

|/jv(x)|i/4 VTV 



dx 
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we have 
(4.69) 



a;- 1 



x + 1 



1/4, , 



|/^(x)|i/4|Ai(-/^(x))|0(-) 



dx 



X sup 
< const C 



-1+5 

cp4>p{cpy + dp) dy 

X(x) 
,{cNCpf'^ 



-1+5 



x + 1 (Arp)i/2 ^ ' 



Also by the analog of the second inequality in H4.fi 5|l for |x + 1| < 5, 



(4.70) 



^i^|Ai'(-/Ar(:E))| <C' 



and hence 
(4.71) 

1/2 



-TV 



-1+5 



a; + 1 



1 



1/4 1 



|Ai'(-/^(x))|o(l) 



X sup 
< const C 



l/iv(x)|V4 
1+5 

cp(f>p{cpy + dp) dy 

X(x) 
,{CNCP)'/' 



dx 



pl/2 



-0{N- 



dx 



.S |x+l|l/4 (iVP)l/2 



i2E£l^O{N-'/'). 



Now we substitute H4.11(l for (pp{cpy + dp) in (|4.fi7(l . Using an obvious schematic 
notation we note that, in view of H4.69|l . (|4.71|l . we have shown 



J" 



: / {Ai + Ai)dx / {Ai + Ai)dy 

Xn ,P -^{^) 

+ / (Ai + Ai')dT / (Ai-0(P-1) + Ai'-0(P-I))rfy 



(4.72) 



oo. 



To estimate the second integral in (I4.72|l we interchange the order of integration. 
Set 2/0 = X{—\ + 5) S (—1 — 5, —1 + 5) and note that X{x) is a 1-1 function from 
[xat.p,— 1 + (5] onto [—1 — (5, j/o]- We conclude that the second integral in (|4.72f) 
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takes the form 
(4.73) 



Li 


y- 


1 




y + 


1 



1/4, ~ 



|/p(y)|V4Ai(-/p(y))0(l) 



y + 1 



y-1 



1/4 1 



l/p(y)P/^ 



Ai'i-fpiy))0(-)]dy 



P 





X - 1 


Ixn,p \ 


X + l 



1/4 , 

\fN{x)\^'^ki{-fr,{x)) 



x + l 



1/4 1 



+ (-l)^+^(c^.cp)l/2 





y- 


1 


'yo \ 


y + 


1 



X- 1 

1/4, . 



\fN{x)\^^^ 



Ai'(-/jv(a;)) da; 



|/p(y)|V4Ai(-/p(y))0(l 



1/4 1 



l/p(y)|i/^ 



-1+5 



,P 



X - 1 



y-1 



1/4 - 



Ai'(-/p(2/))0(i)]dy 



x + l 



X- 1 



1/4 1 



\fN{x)V'' 



k^{-fN{x))]dx 



Note that it follows from g^Hl, 633, 13331, ESZl, that, uniformly for all 

intervals if C [-1 - (5, -1 + iJ], 



y + 1 



l/4 Ai'(-/^(y)) 

l/jv(2/)|i/4 



K 



y-1 



2/+1 



1/4 ~ 

|/jv(2/)r/^ ki{-fM{y))dv 



lUrnil implies that is of order 



This together with 

Thus (cf. H4.72|l ) we have shown (again schematically) 



T" 



(4.75) 



-1+5 



( Ai + Ai' ) dx 



-1+5 



X{x) 



[ Ai + Ai' ) dy 



Now we consider the terms of the form Ai' x Ai' . After changing variables 



(4.76) v^~fpiy)^~apP^^''iy + l)fpiy), u^-~fN{x) 
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and recalling (|4.25|l . I|4.26(l . I|4.27|l and l|4.28|l . we find that the Ai' x Ai' term equals 

(-l)^+^(cjvCp)i/2 times 



(4.77) 



-1+5 



x + 1 



X- 1 



1/4 1 



1 1 



y + 1 



y-1 



1/4 1 



l/p(2/)P/'* 



Ai'(-/p(y))dy 



B]^{x{u)) Ai'(u) 



(5^ap)5/4 (iVP)5/6 

X / Bp{y{v)) Ai'(w) 



1 



d(-1+5) 



Bn{x{u)) Ai' (u) du 



(5^^5^(iVP)^A(_i+,) 

X (^Bp{y{v{X{x{u))))) Ai{v{X{x{u)))) 

- Bp{y{v{-1 + 6))) Ai{v{~l + 6)) 

1 B'pjyjv)) Ai{v) 

^ apP2/3 C/p(y(«)) 

Qi,i + Qi,2 + (3i,3- 



Note from (|4.61|l that a;jv,p + 1 = -(5 + o(l) as N, P^oo, and hence by H4.25|l 



(4.78) 



u{xN.p) +00, u{xN.p) = 0(A^^/^). 



Now 



1 



1 



a, 2 



(aArap)5/4 (Arp)5/6 

X ^Bn{x{u)) Ai{u) 
1 



Bp(y(?;(-l + ,5)))Ai(w(-l+,5)) 

m(2:jv,p) 



m(-1+<5) 
U(.„,P) S^(:c(u)) Ai(u) 



(5;ArAf2/3 7^(_^_^^^ C/jv(a;(u)) 



< 



const 

(iVP)5/6 



since the boundary term and the integrand are uniformly bounded and also the 
length of the interval of integration is 0{N'^/^). Hence 



(4.79) 
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In Qi,3 we integrate Ai'(M) and differentiate the other factors to obtain 
(4.80) 



ru(xjv.p) 



vix(Mu))) B'p{y{v)) Ai{v) ^ \ "^""'"^ 



«(-l+<5) 
u(xn,p) 

u(-l+S) 



-aNN^'^UN{x{u)) 



v(~i+S) Up{y{v)) 

B'p{y{v)) kijv 
v(-i+!>) tjp{y{v)) 



■ dv 



Ai{u) Bn{x{u)) 



d 



X — \v(X(xiu)))\ du 
du '■ ■' 



B'p(X(x{u)))Ai{v{X{x{u)))) 
Up{X{x{u))) 

= Ql,^ + Ql,3 + Ql.S 



where we have used 

(4.81) yov^id. 
Set 

(4.82) l3{u) = v{X{x{u))). 

Lemma 4.5. Uniformly for u(— 1 + 5) < u < u(— 1 — 5), we have as N,P~i-oo 

(4.83) /3(w) = 0(p-i/^) + (l + 0(p-i))M 

(4.84) + 

Proo/. We have in view of 1)4.2611 . (|4.53l) 

-^/3(w) = ^;'(X(x(u))) X'{x{u)) x'{u) 
du 

^ ap CN Up{X{x{u))) 

7V2/3 5^ CP ;7jv(x(u)) 

_ / ^ ^p(X(a;(M))) - Up{x{u)) 

Af2/3 CP yUN{x{u)) Un{x{u)) 

By Lemma El - 1 = 0(P"^), and also X(x(m)) - x{u) = (f^ - l)x(ii) + 
(c/at — dp)/cp ~ 0{P^^) since x{u) G [—1 — (5, — 1 + (5] is uniformly bounded. 
Also by H4.12|l max^gj^j.^+^-j Up{u) < const. In addition, ctN also has a 

complete expansion in powers of 7V~i/(2™) (^gee ()4.12(l 'l. and therefore by the same 
argument as in the proof of Lemma f4. 41 ^ = 1 + 0(P^^), as N,P^oo. Finally 
as Un{x) has a complete expansion (see H4.12II ) in powers of TV"^/'^™)^ uniformly 
for -1 - (5 < a; < -1 + ^, we again find = 1 + O(P-i). Thus 

^P{u) = (1 + O(P-I)) (1 + O(P-I)) (1 + O(P-I)) (1 + O(P-I)) 

= i + o(p-i) 

uniformly for u(— 1 + 5) <u < u(— 1 — (5) which proves 1)4.84)1 . 
Finally by (jI!M)l 

/3(u) = ^(0) + / p'it) dt = /3(0) +u-{l + 0{P-^)) 
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which together with 

/3(0) = viXi-1)) = v{- cn/cp + (dN - dp)/cp) =v{-l + 0{p-^)) 
= -5pp2/3(i _ 1 + O(p-i)) /p( - 1 + 0{P-^)) = 0(p-i/3) 

proves (1133 . □ 

Now we use the uniform boundedness of the functions Ai, B^v, B'j^j, ^j-, and also 
(14.8411 . to conclude the following. Note that the second estimate in ()4.33|l holds on 
the whole of R and that the lengths of the intervals of integration in (|4.80ll are of 
order 0{N^^^). Hence the integral in 3 (and also the whole boundary term) is 
bounded by 0((7V2/3)3/4^ ^ 0{N^^^). The term Q'/g is estimated in a similar way. 
Finally, in Q'/'g we just estimate the integrand by a constant (note (14.841) ). We then 
obtain 

(4-85) < (N^'y (N^"f" = 0(P-') 

I const ./a _ ^ (0/p-2/3N 

l^l,3l ^ p7/3 - (iVP)l/2 

Finally we consider the term in (|4.77() . Using (|4.81() . (I4.82|l we rewrite Qi.i as 
1 1 



(a7v5p)5/4 (Arp)5/6 



«(-l+<5) 



' Bn{x{u)) Ai'(u) 



X Bp{X{x{u))) Ai{P{u)) du 

1,1 



Now in Q'/i, u{xn,p) +00, and we use the estimate for Ai' in (|4.33() (which 
clearly also holds on [—1, +00)) together with the boundedness of the other factors 
to find 

Note that by (j4.83f) . for P and N large enough, 

(4.87) < < C and /3(u) < -1/2 

uniformly for u{—l + S) < u < —1. Hence using (|4.33|l . and recalling that Bn, B'j^ 

are uniformly bounded, and also using the properties of /at after (|4.12|) (cf. (I4.41|l '). 
we find 



(iVP)5/6 7_,-^./3;„(_i+,) 



du < const 



(iVP)l/2- 



We will see below that it is important that we can choose S to be arbitrarily small. 
Collecting the estimates (|OH)l . (jlSHl, (1123, in l|lT7jl and recaUing (|^75)l 
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we conclude (schematically) 

p-l+5 n-l+S 

J44 = / A\ dx Ai dy 

Jxm.p J X{x) 

(4.89) + / Ai da; / Ai' dy + Ai' dx / Ai dy 

XN.P J X{x) ^ XN ,P J X{x) 



Now we consider the term of the form Ai x Ai' in H4.89|l . After changing variables 
as in (|HT)jl and using (glHll), we find that the Ai x Ai' term equals (-l)^+^(cArcp)^/^ 
times 



(4.90) 



-1+.5 



x~l 



x + 1 



\fNix)\'^' Ai(-/jv(x))dx 



-l+<5 



X{x) 



y + 1 



1/4 1 



1 



y-1 

u(xn,p) 



I/P(2/)P/^ 



Ai'(-/p(y))dy 



L]si{x{u)) Ai{u) du 



ti(-l+<5) 
viX{x{u))) 

X / Bp{y{v)) Ai{v)dv 

v(-l+S) 

^ PU{XN,P) 

Ljv(x(u)) Ai{u) du 



aTaT ^'"P"' ^"(-1+^) 



X [Bp{X{x{u)))M{v{X{x{u)))) 

Bp{-l + 5) Ai{v{-l + 5)) 

MX(x(u))) B'p{y{v)) Ai{v) 



dv 



— Q2S + Q2,2 + Q2,3- 



In Q2,3 we split both integrals into the parts with positive and negative arguments. 
Using and iflT^ we find 



const 1 / f du 

\Q2,3\ <T7T7T7^7^ const 



^l/2p5/6 P2/3 V (1 + 1^,1)1/4 



X const + 



(4.91) ^ (l + lwl)!/-* 



" dv 
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Next by 
(4.92) 



|Q2,2| < ^u2p,/6 / ^^(^(") AH") du Bp{-1 + 5) Ai(«(-1 + 5)) 

- (jVP)2/3 ('^^^^^ + ^^^^ Ln{x{u) Ai{u) du^ 0(1) 



r ~1 



(iVP)2/3 



0(1) 



u(-l+S) 



const + ( / Ai(t)dt j Ln{x{u)) 

1 / r' ^ , \ i'Ar(a;(")) , ■ 



since the length of the interval of integration is 0{N'^^^) (again by 14.12|) ) and since 
J'_^ Ki{t)dt and all other functions are uniformly bounded. 

Finally sphtting the interval as in (|4.91l) . and using l|4.12ll . (|4.87() and H4.33|l . we 
find 

|Q2al< '/Tp5/6 ^""^^+ / ' - |Ai(^)||Ai(/3(z.))|dz.) 
const 1 ( 

(4.93) - {NP) 1/2 pi/3 1^°'""* + 7^ ^1/2 J 

(]^7^(--* + (^^''')^'') 
(]^^(^^^^ + ^-^^^)- 



< 



< 



Note that the integral of the form Ai' x Ai in H4.89|l can be estimated similarly after 
changing the order of integration as in (|4.72|l above: the estimate is then the same 
as for the integral of the form Ai x Ai'. Thus collecting the estimates ()4.93|l . H4.92II . 
(|4.91|l in (|4.90(l . and recalling H4.89(l . we conclude (schematically) 



(4.94) 



J44 = / Ai dx Ai dy 

Xn ,P ^i^) 



Note for future reference that we will eventually take a limit S 0, and hence we 
leave only the term 0{S^^'^) in the above estimate (and in similar ones below). 

It remains to analyze the Ai x Ai integral in 14.94|l . As we will see, this is the only 
term that contributes to leading order in J44. Making the same change of variables 
that led to (|4.9n|l . and recaUing H4.a0|l . H4.82|l . we write the Ai x Ai integral as 
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(— l)^+-^^cjvcp times 



{&Napr/HNPy/'\Ju(-i+s) Jo 



Ln(x{u)) Ai{u)du 



(4.95) 



/3(u) 



v{-l+S) 
1 



Lp{y{v)) Ai{v) dv 



Y7^( Qa,! + Q3,2 



(aAr(5p)3/4 (iVP)l/2 



Note that in Q3,2, u > (recall liiiiAr p^^o u{xn.p) = +00). It follows from 
that there is a number uq > 0, uq = 0{P^^/^) such that j3{u) > for u > uq. We 
write 



(4.96) 



+ 



Ln{x{u)) Ai(u) du 

u{xn.p) 



P(u) 



Lp{y{v)) Ai(u) dv 



v(-l+S) 

Ln{x{u)) A\(u) du 



+ 

v(-l+S) Jo 



/3(n) 



Lp{y{v)) ki{v) dv 



Now 



(4.97) 



IQ4 



Lm{x(u)) Ai(u) du 
1 



Lp{y{v))( / Ai(t)dt 
Li,(y(«)) 



/3(«) 
f(-l+i5) 



5p-P2/3 ;7p(y(«)) 



/ Ai{t)dtj dv 



< 



const 

pl/3 



since Uq = 0(P ^/'^), all the integrands in the du and dv integrals are uniformly 
bounded, and the length of the inner interval of integration is 0{P^^^) (recall 

To estimate (54,3 recall first from (|4.27|l that y{v) — (—1) 
Also by H4.12|l . \L'p\ is uniformly bounded. Hence 



apP-^/Wpiviv))- 



(4.98) \Lp{y{v)) - Lp(-l)| < const • \y{v) - (-1)| < const • \v\ ■ p-"^'^. 
Thus, since /3(w) > 0, using 14.26|l . H4.33|l and the fact that /Ar(— 1) = 1, we obtain 



/3(«) ./3(u) 

Lp{y{v)) A:i{v)dv ^Lp{-1) / Ai{v) dv + O 



I3{u) ^g-(2/3)^='/^ 



p2/3 



■ dv 



Pin) 



2^/4 / Ai(w)dt. + 0(p-2/3), 
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and so 



(4.99) 



^4,3 



Pin) 



Ln{x{u)) Ai(u) 2i/4 / Ai{v) dv + 0(P-"/3) du 



=2^/4 



u{xm,p) 



/3(u) 



Ljv(a;(u)) Ai(u) ( / Ai(i;) dw) + (^(p-^/a) 



again by H4.33|l . Next 
(4.100) 

u(xn,p) 



L N {x{u)) Ki{u) du I Pl\{v) dv 



I3{u) 







Ln{-1) 

u{xn.p) 



u(xm,p) 



/3(u) 



M{u)du / Ai(i;)dv 



+ / {LNix{u)) ~ Ln{~1)) Ai{u)du Ai{v) dv = Q^^i + 



/3{u) 



15,2- 



Now by the uniform boundedness of and 14.27|l . and since \v^ Ai{v)\dv < oo, 
J =0,1, 



/•u{xn,p) 
IQ5,2| < C / 

J Uo 

Next 
We have by 



«(a;jv,p) 



|x(u)+l| I Ai(u)| du < ci7V-2/3 / |u Ai(u)| du = 0{N'^/^). 



"(a;jv,p) 



Ai(u) du 







Aiiv) dv = Ql, + Q'' 



Ai{v) dv < c\P{u) -u\ = 0{P-^/^ + VI) 



and hence 



/>C30 

|Q'5'i|<2i/4 / |Ai(u)|0(p-i/3+p-i|u|)dM-0(p-i/3). 
Jo 



Also 



%,,^2-(/ 



Ai(u) du / Ai(u) 



Jo J u{xn,p) . 

= 2i/M Ai{u)du Ai{v)dv + 0{p-^/^) + Oie~'^) 
Jo Jo 

where we have used the fact that uq = 0{P~^^'^), the uniform boundedness of 
/p" Ai{v)dv for m > 0, and the (super)exponential decay of Ai(i;) for v > 0. Thus 



(4.101) 
N, P-^oo. 



^4,3 



+ 00 



2^/2 / Ai{u)du / Ai{v) dv + 0{p-^/^), 
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To estimate Qa^2 we first write 



Lp(y(v)) Ai(v) dv 



(4.102) 



= Lp(-l) / Ai{t) dt - Lp{y{v{-1 + S))) 



Ai{t) dt 



(-1+5) 



/ Up{y{v)) 



dv. 



But by as v{-l + 5) = -apP^/^5fp{~l + 5)-^ -oo, 

rv(-l+S) 



Lp{y{v{-l + 5))) 



M{t) dt 



< const \&pP^^^6fp{-l + (5)1-3/4 
1 



< const 



and also 



(4.103) 



Thus 



(-1+5) 



Ai{t) dt 



^3/4pl/2 

L'p{y{v)) 



Up{y{v)) 



dv 



< 



const 



Q4,2 



p2/3 

< const 

I 

Ln{x{u)) Ai{u)du 







dv 



&pPV3Sfp(-i+s) (1 + 
51/4 



pl/2- 



Lp{-1) j^^ A\{v)dv + 0[j^ + 



1 \ 



pl/2 ^3/4pl/2, 

n{x{u)) Ai{u)du [ Ai{v)dv + 0(j—^^^ 



21/4 



where we have again used the estimate in (|4.33|) as in (|4.99(l . Now as in the analysis 
of KTm we find 



Ln{x{u)) Ai{u)du 



tl(xjv,J=) 



Ai(M) du + 



ti(a;jv,p) 



{Ln{x{u)) - Ljv(-l)) Ai(u) du 



= / Ai(M) du + 0{P-^'^ + e-"^) + 0{N-^'^). 
Jo 

We conclude that 

r+oo /-O / 1 \ 

(4.104) g4,2-2i/2y^ A\{u)duj^ Ai(«)dz; + 0(p-V3) + o(^___j. 
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Thus as N,P^oo, the second integral m (|4.95|l behaves as 



(4.105) 



i3.2 



+00 



Ai{u)du / Ai(«)dw + 0(p-i/3) + 0(,5-3/4p-i/2)_ 



The first integral in (|4.95|l is given by (recall l|4.81(l 'l 
(4.106) 

f/3(«) 



Q 



3,1 





«(-l+(5) 



Ln{x{u)) Ai{u) du 



Lp{y{v)) Ai(u) dv 



v{-l+6) 



{x{u)){ Ai{t)dt)( Lp{y{v))ki{v)dv 



1 



Ai{t)dt 





«(-l+<5) 



du 



' \Jn(x{u)) 
X ( / Lp{y{v))Ai{v)dv) 

Ai{t)dt) Ln{x{u))Lp{X{x{u))) Ai(/3(u))/3'(M) du 



' u{ — l+5) J —OQ 

= Qea + Q&..2 + Qe.s- 
Observe that u(— 1 + 5) < u < 0, and so f3{u) < const < cx3 by (|4.83|l . 



/3(«) 



■u(-l+<5) 



Lp{y{v)) Ai(7j) c?w 



Lpiviv)) 



(4.107) 



< 

v{-~l-\-S) J —OQ 



Ai{t)dt 
Ai{t)dt 



I3{u) 
v(-l+S) 

L'p{y{v)) 



Up{y{v)) 



■ dv 



< const + 0((iV2/3)i/4p-2/3^ < ^Qj^g^ 
where we have used the bound on the integrand (|4.34|) . Thus again using (|4.34|l . 



(4.108) 



const 



C 



7V2/3 (1 + |u|)3/4 



du = 0(7V-i/2). 



Note next that in 



Qe,! =ijv(a;(0))(y Ai(i)dt 
-Ln{x{u{-1 + 5))) 



m 

Lp{y{v)) Ai{v) dv 

vi-l+5) 

ui-l+S) .piui-l+S)) 

Ai{t)dt) ■ 



Lp{y{v)) Ai(u) dv 



the second term is 0{P-^/'^) by H4.107|) . As in H4.102|l . the first term 

can be written as 



Aiit)dt -Lpiyim)) 



0) . 

Ai{v)dvj +o(r3/4p-i/2) 

= 2^/2 ( r Ai{t)dt^ ( r Ai(,;) di;) + 0(<5-3/4p-i/2 ^ p-i/3) 
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in view of fOTjl and Thus 

(4.109) Qg^i = 2^/2 Ai{t)d?j +0((5-3/4p-i/2^p-i/3)^ 

Finally, 
(4.110) 
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■a 



m(-1+<5) 



Ai{t)dt^ Ln{x{u))Lp{X{x{u))) Ai{(3{u)) P' (u) du 

Ai{t)dt] Ln{x{u))Lp{X{x(u))) Ai(/3(u)) du + 0{p-^'^) 



by (|4.84() . the uniform boundedness of the remaining factors in the integrand, and 
since the length of the interval of integration is 0{N'^^'^). Denote 



A{u)= [ Ai(/3(s))( / Ai{t)dt]ds. 



Si 



mce 



du 



A{u) 



Ai{t)dt] Ai(/3(u)) 



the remaining integral in 14.11()|l after integration by parts becomes 



-A{u) Ln{x{u))Lp{X{x{u))) 



u(-l+S) 



(4.111) 



— A{u) ^^4^ Lp(X{x{u))) du 
unN^I-^ Ju(-i+s) Un{x{u)) 

1 A( \ T ( ( ,, L'p{X{x{u)))X'{x{u)) ^ 

~ ..2/3 / A{u)Ln[x{u)) , , du. 

aNN^/-^ Ju(-i+s) Up{x{u)) 



We need the following result. 

Lemma 4.6. As N^P^oo, we have uniformly for u(— 1 + 5) <u <{) 

(4.112) A{u)^( Ai(s)( / Ai{t)dt\ds + 0{5^''^) + 0{5-'^''^N~^). 

J — oo J — oo 

Proof. For m(— 1 + 5) < u < — 1 we recall (|4.76|l 

(4.113) u{-l + 5) = -aN5fN{-l + 5) = -c5N'^/^{l + o{l)), c> 0. 
For u < 



(4.114) 

Note that 
(4.115) 



Ai(/3(s))f / Ai[t)dt) ds 

,(-l+<5) ^J-oo ^ 

Ai(s)( / Ai{t)dt 



u(-l+S) 



Ai(s)( / Ai{t)dt]ds 
ds 



"A 



u(-l+<5) 



i(/3(s)) - Ai(s)] ( / Ai(t)dt) 



y Ai(s) ( J Ai(t)dt) ds^^(^J Ai(s) ds 
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and so, by I|4.1ia|l . the first term in the RHS in H4. 114(1 equals 

i( / Ai(s)ds) - o((57v2/3)-3/2) = o(r3/2iv-i) 

In view of H4.22|l , (|4.83|) , we bound the second term in (|4.114|l by 

ds 



«(-l+<5) 



Ai(/3(s)) - Ai(s)| • / K\{t)dt 

J — OO 

( sup |Ai'(i)|) •|/3(s)-(s)|. K^(t)dt 

«(-l+(5) ^tG[/3(s),s] ^ J-oo 

1 



< 



< const/ (l + |s|)i/4 0(p-i/3) + |s|0(p-i) '^^ 



ti(-l+(5) 



l + |s|3/4 



uniformly for u(— 1 + (5) < w < 0. □ 



From H4.115|) . Lemma Ol (|4.111|l . and the fact that A(u{-\ + (5)) = 0, we 
conclude that 

(4.116) Qg3^__M Ai(s)ds) +0(<5i/2)+0(r3/2iV-i) 

again by the uniform boundedness of the terms in the integrands, the integrability 
of (/"^ Ai(s)ds)^ - u ~* -OO, and the fact that u{-l + S) = 0{N^/^). 

By bA . imni, (ICTKll we see that KTm is given by 



(4.117) 



Q3,i = / Ai(s)( / Aiit)dt\ ds 

+ 0(<5l/2) + 0(,5-3/47V-l/2) + 0(5-3/2^-1) ^ 0(p-l/3)_ 



Thus by 623, l|4.117|l . (|ni751) . the Ai x Ai integral in equals (-1)^+^ {^^t^Y/'H^Npyn 

times 

(4.118) 

Ai(.) ( m)d^ ds + 0{P-^'^) + O (5V2 + _L_ + 
1 /•+00 2 / 1 1 \ 
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by glEl), and since Ai{t)dt = 1 (see |XbSt] ). Thus from ^J^{2), 



J44 — J44 + J44 — 1) 
1 



N+P 



Am 



C)(p-l/3) +0(p-l/(2rn)) 

+ 0((5i/2) + 0{P-^/^) by El??!! 



by 



1 + 0(P-1/3)^(9(P- 



l/(2m 



4m (iVP)i/2 

+ 0(51/2 ^ ^-3/4p-l/2 ^ J-3/2p-l^ 

By this estimate and (|4.57|) . I4.58|l . (|4.59|l . we find finally that 



(4.119) 



(-1) 



JV+P 



(-1) 



N 



4m 



2™ 



0(^1/2 +p-l/2^-3/4^p-l^-3/2 



In a similar way we obtain 

^1+5 

J2 = 



.1+5 .1+5 

' cn4'n{cnx + cIn) [ / cp(j)p{cpx + dp)jdx 
1-5 ^Jx(x) ' 



(4.120) 



(c^cp)i/2 
(Arp)i/2 



— +0(P-l/3)^0(p-l/(2m)) 

4m 

+ 0(51/2 ^ p-l/2^-3/4 ^ p-lj-3/2^ 



(Note that the integral leading to the term in (|4.119|l . is not present for J2.) 

It remains to consider J3. Again assume that N,P~^ 00. By 



/• r+oc 

■h^CNCp j (j)N{cNX + dN)dx / (j)p{cpy + dp)dy 
Jh Jx(x) 



1-5 



CNCp / (l)NicNX + dN)dx[ / + / + / ]<t>p{cpy + dp)dy 

JI3 \Jx(x) JI2 Jh 

J33 + "^2 + </31- 



By (I33HI, iFT7|) . 

|J32l^l7„.l_.l/2 



I + |J3,| = 0{{cp/PY'^ + {cp/Py/^ e-^-) = 0(A^-V2). 



To analyze J33 we first introduce p G /a such that 

X{x% p) = x% pCN/cp + (rfjv - dp)/cp = ±(1 - S) 
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which in view of (|4.54|l . imphes that for N,P^(X), 

(4.121) X{x)el3, xel'3 = [x]^p,x+p]cil3. 

Set /g = /a \ and note that by Lemma lOI 



(4.122) 
Then 



p = ±(l-5) + 0(p-i). 



J33=CNCpyJ +1 ](j)N{cNX + dN)dx 
= V + J" 



1-S 



X(x} 



(jjpicpy + dp)dy 



Now 



(4.123) |4'3| < / ^ 



1/2 

Cat 



1/2 ^1/2 '^^ - 1^3 I - (-^p)l/2 ^ 



by 623, Bl5ll . (TOa . 

Now taking into account H4.121|l and using H4.13|) 



4 



33 



-y/CNCp 



dx 



1-S 



|1 -a;2|i/4 
dy 



(cos/+(a;)) (l + 0(7V-i) + (sin/^(x))0(iV-i) 



where 
(4.124) 



(cos/+(2/)) (1 + 0(P-1) + (sin/p(y)) O(P-I) 



TV 



|1 - y^|^/^/iAr(?/)dy ± - arcsincc 



and /ijv satisfies (|4.7f) and H2.7|l . 

Clearly a product of two terms, one of which has a factor 0{N~^) and the other 
0{P~^), gives rise to a smaller order contribution ^^j^^ 0{N^^). Consider next 
the integral 



dx 



|1 -a;2|i/4 

< 0{N-^) 

< 0{N-^) 
= 0((iVF)-i), 



(sin/-(x))0(iV-i) 



1-5 



■; 777TTT I IllcLX 

|l-x2|l/4y _l+5<:ri<l-d- 



dy 



|l-y2|l/4 



(cos /+(?/) 



where we have used (|4.21ll in the last step. Thus 



-y/CNCp 



dx 



|1 -a;2|i/4 

_ ^CNCp 



'NP 



sin/^(x))0(7V-i) 



0(iV~i). 



1-5 



dy 



(cos/^(?/)j 
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Similarly, interchanging the order of integration as in l|4.7l'i|l , we find 



^\fcNCp I cosf^{x) 



dx 



1-:e2|i/4 



l-<5 



(x) 



sin/p(y))0(p-i)- 



dy 



-\JCNCP 
TT 



y2|l/4 

(sin/p(y))0(p-i) ^ f^-'(y^ cos/+(x) 



1+.5 



|l-y2|l/4 



■ dy 



l+d 



|l-a;2|i/4 



dx 



x/NP 



0{N-^] 



Thus 



(4.125) Jg = -75^ 



cos/;^(a;) 



da; 



1-5 



Note that, in view of (|4.7() . for a fixed 5 > 0, there exists 775 > such that as 

N (X) 



(4.126) 



(/^(2;))' = ^|l-:^'l^/'/iiv(x)± 



1 



uniformly for x E I'^. The integral in (|4.125|l is given by ^^/cnCp times 



cos/j^(x) 



sin/+(y) 



(/^(y))'|i-y^p/^ 



l-<5 



X{x) 



(4.127) 



l-<5 



1 



\i-y'\'/'ifMy)y 



dy 



(sin/+(y)) 

sin/p(-^(^)) da; + 0((iVP)-i) 

ii - -'\'/' frix{x))\i (x(x))2|i/4 + ^^^'"^^ ^ 



in view of (|4.12t)|) , since we can integrate by parts in x once more for the boundary 
term at y = 1 — 5, and also integrate by parts again in the integral term Jx(x) ' ' ' ■ 
Note next that by Lemma 14.41 we have again uniformly for x E 



\X{x)-x\=0{p-'). 

Also, uniformly for x £ I3, \fp"{x)\ < const • P and |1 
X{I^) C I3, we conclude that 

1 1 



< c. Thus, as 



f+\Xix))\l - (X(x))2|l/4 /+'(a;)|l - x2|l/4 



< 



const 



and hence using (|4.127|l we find 



(4.128) J3 = - -Vq^ 



cos/+(x) sin/+(X(x)) ^/cnCp 



n /+'(x)|l-x2|i/2 



da;+ 0(jV~^/^) 



Recalling (|4.12t)|l . (|¥7|) and (|^ . we derive from (|4.128|) 
(4.129) 



= 



4 ycjycp I" cos/+(a;) sin/+(X(x)) ^ ycjycp ^^^r-i/(2m) 



(1 - a;2) /^(a;) 



'NP 



0{N- 
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Now the integral above is given by 
(4.130) 

cosf+{x)smf+{X{x)) , 1 f sm[f+iXix))+f+ix)] 



{l-x^)h{x) 



■ dx — — 
2 



dx 



{l-x^)h{x) 
1 f sin[/+(X(x))-/+(x)] 



il-x^)h{x) dx^Qr,^ + Q,,,. 

Integrating by parts in Qy^i in the same way as in ()4.127f) and noting that for 
sufficiently large N, by H4.126II . 

f+'{X{x))X'{x) + f+'{x)>Nf,s 

uniformly for x £ I^, we conclude that 

IQt.iI <o(iv-^). 

Finally consider Qt,2- Introduce the convenient notation 

1 



(4.131) 



6ix) 



\1 - t^\^^^ h{t) dt 



and note that 9{1) = tt/2 by (|2.5(l . H2.7|l . the fact that h{t) is even (see H2.8|l l and 
the normalization condition J^^ 'il}'^^'^\x) dx ^ I. 

Lemma 4.7. Uniformly for x £ as N, P^oo, 

0(7V-i/(2™)). 



(4.132) /+(X(x))-/+(a:) = (iV-P) 
Proof. We have 



±^xO'ix)~Oix) + l 



f+{X{x)) - f+{x) = [fUX{x)) - /+(a:)] + [/+(x) - /+( 



F1+F2. 



By mm 

F2 = \f\l- v^\'/^[Php{x) - NhM{y)]dy. 

Note that by ifTTjl with q = 2m and (|^?5T|l 

Php{x) - Nh^ix) ^ P{hp{x) - hN{x)) + {P- N)hN{x) 

2m 
k=l 

+ (P - N)[h{x) + 0(7V-i/(2'"))] 
= (P-7V)/i(2;) + 0(iV-i/(2™)) 
uniformly for —1 < a; < 1. Hence, uniformly for x £ I'^, 

F2 HP - N)^ |i " y^\"My)dy + o{N'^'(^"^'>) 

= -{N - P) [e{x) y^\"^h{y)dy\ + 0(7V-i/(2™)) 

= -(iV - P) [6l(a;) - |] + 0(iV^^/^^"^), 
since 9{l) = 7r/2. Next 

Pi =f+'{x) {X{x) -x)+ fnm) {X{x) xf/2 
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for some ^(x) between x and X{x). Note that by (j4.53|l and Lemma lOI 

X{x) -x = {cn/cp - l)x + (djv - dp)/cp - ^^^^x + 0(^-1-1/(2")). 

Also by H4.126|l and \fp"{x)\ < const • P uniformly for x e I3. But 

P, 
"2 



/^'(x) = -|1 - x2|l/2/t(.T) [1 + C)(p-l/(2™)) + O(F-l) 



and hence by (|4.131|) 

2 2™ 2m 
which completes the proof of the lemma. □ 

We need the following fortunate and remarkable fact. 

Lemma 4.8. For any V(x) — K2mx'^"^ + - ■ ■ the function 9 defined by H4.1c{l|l solves 
the linear first order ODE 

Oix) x9' (x) = arcsinx, —1 < x < 1. 

2 m 

Proof. See Remark l6.1l in Subsection l6.3l □ 
Thus from (|4.130|l . (|4.132|) 

J 2 r sin [{P ~ N) arcsinx + {N - P)f )] ^ ^,^l/(2m)^ 

^ TT y/NP Ji;^ (l-x2)/i(x) ViVF 

Note that by (|4.122|) we can now replace /g with = [— 1 + (5, 1 — (5] introducing 
an error ^^yj^ 0{N~'^). Now using the formula for the sine of a sum and noting 

,1 J sin[(P— A^) arcsin a;l • ij r n :i 

that — 4n — 'ii 1 r \ — - IS an odd function we hnd 

(1 — 2;^) h{x} 

2 ^/qTcF ^^.^^ / (A^ - P)Tr\ f^^^ cos [(P - N) arcsinx] 



TT \/]VF V 2 / (1 - x2) h{x) 



Assume iV - P is even. Then J3 = 0(iV-i/(2™)). Assume TV - P is odd. 



Then 

cos [(P — iV) arcsinx] 



O , x-±(l-0) 



(l-x2)/l(x) V|l -x2|l/2y' 

and hence for N ~ P odd 

2 Vcjvcp /(A^-P)7r\ /"i cos [(Af- P) arcsinx] , 

J3 = . sm j / . . ^ ^cia:: 



(4^^33) TT VAP V 2 7 7-1 (1-.t2)/i(x) 



^/]VP V ' 
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By (jl3Zj), (11321), gSni), as N,P^oo, 

(4.134) 

(0Ar, e(/)p) 2 ( / dxj / (/)p(y) dy) 



CNCp( + )(f>N{cNX + dN)dx / (f>p{cpy + dp)dy 
^Jli Jl.,'' Jx(x) 



Ji + J2]+Ji 



^^/i ^ (1 + + 0(iV-l/2) + 0(^-l/(2H)) 



X (1 + (-1)^ + 0(p-l/2) ^ Q(p-l/(2m)-)) 

0(e-"^^) by 



(-1) 



(-1) 



AT 



1 



4™ 2™ 4m 



(7V-P)7r\ 2 /"i COS [(A^ - F) arcsina;] 



( 2^J 



(l-a;2)/i(x) 



dx 



\/CNCp 

Vnp 



+ 0(iV-l/(2m)) ^_ (9(^1/2) 

(-l)^-(-ir 



by 



4m 



+ /(AT _ P) 



+ 0{P-'^') + 0(p-l/(2™)) + 0(^jl/2 ^ p-1/2^-3/4 ^ p-lj-3/2j 



where 



s(g arcsin a:;) 



dx, q odd 



^i"^ ('2 ) ^ /-I (i-K^) ^(k) 
0, q even. 



As {(j)N,^4'p) = — (e0A, '/'p) = ~(0P, c'/'a) we see that (|4.134f) is true also if < P. 
Hence as N, P^oo 



Vc^i^^-I{N-P) + 0{5^^^) + os{l), N-P odd 
^/NP [0(51/2) + 05(1), TV -P even, 



where 05(1) = 0{P-^/^) + 0(P-i/(2™)) + 0(P-i/25-3/4 + p-15-3/2). Note that 
some of the smaher order terms in 05(1) are proportional to for some a > 0. 
Nevertheless we see that e.g. for — P odd 



-0((5i/2) < liminf 

N,P^oo 



NP /(-l)^ 



< limsup 

N,P^oci 



y/CNCp 

Wp 



^CNCp 



2m 



{e<j>M,<l^p)-[^-^~I{N^P) 



< 0(5l/2). 
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Letting (5 — > we conclude that lim inf = lim sup and hence 

■(-1)^ 



(4.135) (e</.^,0p)^^/^ 



2m 



/(iV-P) + o(l) 



VNP 

as N, P^oo, N ~ P odd. Similarly 

(4.136) (e^^,0p)^^po(l) 

as N, P^oo, N — P even. This completes the proof of Theorem 12 .21 

Remark 4.3. As noted above, the proof of Theorem 12.21 involves a lim inf /lim sup 
argument as N,P^oo, followed by a limit as (5 — > 0. However, as noted at the 
beginning of Subsection 14.2.11 the error estimates from |DKMVZ2] in (|4.6|l et seq., 
are uniform only for S in compact subsets of (0, Sq] , and the reader may be concerned 
that the 0(6^^'^) term in the lim inf /lim sup argument above, in fact depends on 
constants, say, that blow up as ^ — > 0. But the reader may easily check that all 
the error terms that arise in the proof of Theorem 12. 21 and that could blow up as 
5—^0, are always multiplied by negative power of N or P, and hence vanish in the 
lim inf /lim sup argument: all that remains are 

0(5^/2) tg^j^s that arise from the 
evaluation of integrals with explicit (5-independent integrands (see e.g. the proof of 
(I4.133|l above), and do not blow up as (5 —> 0. 

Remark 4.4. In the case m — 1 the OP's are just the Hermite polynomials. They 
have the exceptional property that if N is odd then ecfy^ is a polynomial of degree 
— 1 times and hence (e^jv, 'i>p) is identically zero for N odd and N < P. 

However, we have 

cos (q arcsin x) [tt, q— 1,5,9,- ■■ 

dx — { 

l"a;2 |-7r, g = 3, 7, • • • , 

and for m — \, h{x) = Ahy H2.8|l . A simple calculation now shows that the leading 
coefficient in (14.13511 is zero for A^ odd, N < P, P even, so that the calculations 
match. 

Remark 4.5. The following (nonrigorous) argument is consistent with our asymp- 
totic formulae. We know by Theorem 12.11 that for any polynomial V the matrix 
Doo is banded and looks asymptotically like the product of a Toeplitz matrix, say 
Doc, whose diagonals are given by certain binomial coefficients times the diagonal 
matrix Too = diag((mK2m&j'"~^)i>o)j bo = 1 (see H2.2|) '). Take a large enough even 
N and let Dff he an N x N section of Doo, i.e. {DN)i,j = iDoo)i,j, < z, j < A^— 1. 
Also set 

Tn = diag((mK2™6,'"-^)fj-o\ &o = l- 

Consider a submatrix in the middle of the matrix {DnTn)~^ near its diagonal 
which is small compared to A^. For large A^, in view of Theorem 12. 5f il. we would 
expect this small submatrix to look like the corresponding submatrix of the matrix 

Let us take, e.g., m = 2 and V{x) = x^. Then n = 2m —1 = 3. In this case the 
Toeplitz matrix Dn has diagonals (cf. Theorem 12.111 

(••• 1 3 - 3 - 1 •••)■ 
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In view of the preceding discussion, we would expect that for large enough even N 
the submatrix of T^^D^^ formed by the rows M — n, - ■ ■ , M — 1 and the columns 
M, • • • , M + n — 1, Af even and M ^ N/2, to look like the submatrix of Eoq formed 
by the same rows and columns. In other words, we would expect that this small 
submatrix of would look like the small submatrix of Tj^eN located at the same 
position. (Recall that the rows and columns in DJj^ and Tn^n are enumerated 
from to iV — 1 so that the submatrix chosen above corresponds to B12 defined in 
(fO^ about the center (M, M).) 

The overall coefficient multiplying the elements in the middle of the matrix T/vCat 
is the same as ^J^, that is mK2mb'^M~'^ cmM'^ = ^^^(1 + o(l)), for M large 



enough. Using Maple software we have computed the matrix Dj^^ for = 20. It 



turns out that the submatrix of Dpj with rows 7, 8, 9 and columns 10, 11, 12 (here 
M — N/2 = 10, n = 3) and the corresponding submatrix of TArejv computed using 
Theorem 10 fhere h{x) = ^{1 + 2x^) hy are given by 

/-0.01630 0.00435 \ /-0.01635 0.00438 \ 

0.15078 ' 0.15032 

V 0.06113 -0.01630/ V 0.06100 -0.01635/ 

respectively. We see that the corresponding nonzero elements differ by at most 
0.8%. For N = iO, M = N/2, the maximal difference is aheady less than 0.002%. 
We have also done similar computations for m = 3, 4, n = 2m— 1 again for N — 40, 
M — N/2. It turns out that the nonzero elements in the small blocks computed 
using the matrix inversion and Theorem 12.21 as described above differ by at most 
0.005%, 0.02% for m = 3,4, respectively. 

4.3. Convergence of derivatives and integrals of the Christoffel Darboux 
kernel for weights e~^^^^, V{x) = K2mX^"^ + '^2m > 0. We start with 
the convergence of the derivatives in the 12 entries of the kernels Kn,i3, (3 — 1,A. 
Our main results in this direction are (|4. 150(1 . ((4.151(1 and their Corollaries ((4.157(1 . 
(14.158(1 below. After that we prove the convergence of the integrals in the 21 entries 
of the kernels Km,p^ P = 1,4. The main result for the integrals can be found in 
(14.171(1 below. 

4.3.1. Derivatives. Fix Lq > 0. For r, ^, 77 in the compact set |^|, |r| < 2Lo, 
define 

(4.137) Dr,N{^,v) = —KN(r+^,r+^ 

QN \ QN qN 

In what follows, qn is any sequence of numbers with the property 

(4.138) iN^^^^il + OiN-'^)) 
for some < a < 1, as iV — > cx). 

Remark 4.6. In the calculations that follow we use formulae for Z3r,jv(^7?y) based 
on ((1.7() that holds for ^ ^ t]. But Dr,N{£.,v) is continuous in ^,7/ (see e.g. ((1.511 1. 
and it will be clear from the calculations that we can obtain analogous formulae for 
Dr,N on the diagonal by taking the limit 77 — > ^ in the formulae below. 

Eventually in 1(4.1371) . ((4.138(1 we wiU take qn = RN,is{r),RN/2,iAi'^) fo'^ P = 
1,4 respectively, but at this stage asymptotics for i?Af,i,i(?'), i?Af/2,i,4(?') of type 
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KT^ must Still be prove d, see ^4.152|l . ^JSM, below (of. D, p. 240 et 

seq.] where the analog of H4.138|l is proved for Rn,i,2{0) — ^^^(0,0) in the case 
V{x) = x^™ without lower order terms). By p. 3(1 . qjq ^ j^^-'^li'^'^) to leading 
order. 

We use the notation in |DKM V!Z2] . In particular set 



g(z) = gjv(^)- y \o^{z~x)d^^'''^\x)^ j ^\og{z ~ x) - x""]^^^ hN{x) dx, 
z e C \ (— oo, —1], and for z e (—1, 1) let 



-AT 







S+iz) = I I e--^''Y+icNZ + dN) 

(4.139) 








^ ^-Nig+{z)-l)as 



here ± refer to the boundary values from above/below the real axis, respectively 
(cf. |DKMVZ2[ (4.22)]). Here Y solves the Fokas-Its-Kitaev Riemann-Hilbert 
problem fot the polynomials orthogonal with respect to the weight e~^^^^dx (see 
pKMV Z2 Thm. 3.1], and the constant I = In is given by (5.35), loc. cit.). Finally, 
for z € (-1, 1) let 

(4.140) ^n{z) = g+{z) - .g_(z) = t f |1 - x^l^/^h^ix) dx 

J z 

and set 

rN ^ {r - dN)/cN- 

Note that ^jv(z) e and tn = 0(7^1/(2 ")) by fT^ . fTH) . We need the 

following properties of 5*+ which are proved in |DKMVZ2l Sec. 7]: 

(i) For -1 < z < 1, 

(4.141) detS'+(z) = l. 

(ii) Let a+{z) ~ (fijrY)^* and set 



^ I a+(z)+a+(z) ^ i(a+(z) ^ - a+(z)) 
i(a+(z) - a+(z)-i) a+(z) + a+(z)-i 



^f'(z)^-| I, -l<z<l. 



Then as N ^ oo, S'+(z) converges with all its derivatives to 5^''(z) for all z in 
compact subsets of (—1, 1). In particular, for any < (5 < 1, 



fc 



d 

(4.142) sup max — ^S'+(z) 

for fc = 0,1,2, - • •. 



< Cfc < oo 
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Now a simple calculation (cf. (Dj p. 24]) shows that 
(4.143) 

= ( 5+11 frAT + ^^'je^^" +5+i2frjv + — ^')e-*«"('^"+''/(«""")) 

X ( S+2i(rN + e^«"('-"+«/(?"^")) +S+22(rN + ) e'^^" 

where ^ rj) indicates the same terms with ^ and rj interchanged. Consider the 

terms 

(4.144) 

V QnCN'' ^ qNCN' 

S+iiirN + ^—1 5+22 frjv + 1 («"(""+«/^''"""»-^"(""+''/(«~""») 

V qNCN^ V g^CA^' 



5, 



/(2^i(C - ry)) = gs.i + Q8,2 + Q8,3 
where 

5+ii(rjvH — ) — 5+ii(rArH ^ — ) 

„ _ ^ V qNCN ) \ qNCN ) 

X 5+22 f?'A + — ) e'T{iN{rN+v/iqNCN))~iN{rN+i/{qNCN))) 

and 

Qg 2 =5+11 (rN + ^ ) g-f(?N(l-N+»?/(?NCjv))-Cjv(rN+?/(9JVCN))) 

27r*(e-r;) 

and also 

Qs.3 =S+ii(rN ^ — ) 5+22 fr-w H 

V qNCN'' V qNCN'' 

g-j- (?N(rjv+?;/(givciv)) — CN(rjv+?/(givCN))) _ g"T (r jv +4/ (g]vcjv)) — ^iv(rjv 

2^*(e - ^) 

Note that 
(4.145) 

^+ii(rA + ^) -^+ii(rA + ^) _ 1 /-i^^ / ^ ?? + T(^-ry) x dr 
2Tri{£_-ri) ^ni Jq V g^cjv ^qNCN 

and 

f f^A(rA + ^^)-eA(r-A + ^)) =*^(e-'7)G^(^,'7) 
(4.146) 2 V g^CA gwcw / 
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as — > (X), where we have denoted 



and observed that by H4.138|l . 

(4.148) GNit V) = (1 + 0(A^-^)) = 1 + 0(jV~°), iV ^ oo. 
From the formulae above and the relation '^(^"''^ = , , we see that Qs i, 

2(3ivCiv hjv(0)+o(l)' ~*o,ij 

and similarly (58,2, and all their ^,77 derivatives converge to zero uniformly for 
1^1, |7?|, \r\ < 2Lo, with errors of order 0{N-'^). 
We also have the formula 

Qs,3 = 5+11 frjv H — ) S+22(rN H ^ 

(4.149) ^ '^"'"^ ^ 

xGNitv) cos (T7r(^- 77) Gat (^,77)) dr. 

Together with similar formulae and calculations for the other terms in ()4.143|l we 

see thati?,,jv(e,^) = +Q(A^'^) = Xoo(e, ^7) + 0(jV-") + 0(jV-^) 

where -R'oo(C ~ ^7) is the sine kernel The same is true for all the ^,77 

derivatives. More precisely we have for any j, k > 0, 

(4.150) ^^—KJr+^,r+^]^^^K^{£,^^) 

uniformly for I'CM'yIikl < 2Lo as iV 00. Note that the terms of the form 
5'+ii(- • • )S'+2i(' • • ) and 5+i2(- • • )S'+22(' • • ) do not contribute to leading order, and 
also S+ii{rM)S+22{rN) - S'+i2(7'7v)<S'+2i(fAf) = 1 by (|4. 141(1 . Moreover, keeping 
track of the estimates, we find that 

(4.151) the error term in (|4.150ll is 0(iV"") uniformly for hi, \r\ < 2Lo. 

We see from the above that the largest error in (|4.152|l arises purely from the 
asymptotic evaluation of ^^^^"^^ using l(4.138|l (cf. I(4.148|l ): apart from this term 
the error in H4.152|) is 0{N-^), rather than 0(7V""). 

4.3.2. We can now verify asymptotics for KN(r,r) of type 1)4.138(1 . Indeed from 
the uniform convergence above, taking = for which we may, in turn, 

take a = 1 in ((4.138(1 . we have for \r\ < 2Lo, as iV — > 00, 
(4.152) 

KN{r, r) = qN Dr.N{0, 0) = qN lim Dr,N{L 0) = 9w f lim "'".^ + 0{N-^)) 

C^o \C^o 7r(4 — U) / 

= q^(l + 0{N-^))^^:^{l + OiN-^)). 

This means in particular that ((4.150() . 1(4.151(1 are true for g^r = KN{r,r) with 
a = 1 for any \r\ < 2Lq. Hence 

(4.153) KN{r, r) ~ ivi-i/(2m)^ foj. |^| < 2L0. 
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Estimating the correction terms for Sn. lir, r) as in (|3.5|l . we find for any |r| < 2Lo 

(4.154) ^fLpll^i + o{N-'/'). 

KN{r,r) 



Similarly for any \r\ < 2Lo 
(4.155) 



Knir, r) 

Note that (|4.154l) . I4.155|) together with H4.153|) prove (|1.23f) . and hence p.24|l . 

4.3.3. Derivatives, continued. We now consider the term 

(4.156) XN{i,v) = \{KND){r + -^.r+^], ^at as in glSHI)- 

\ IN qN J 

Terms of this form arise in the 12 entry of ^if^"'^''(r + + Clearly 
XAr(^,7y) = -^Dr^Nii^ri) (cf- H4.137|) ) and so by 14.150(1 we see that H4.156|) satis- 
fies 

(4.157) XN^i.v) - ~5;^^oo(e - ??) = ^^oo(^ - r?) 

together with all its ^,77 derivatives, uniformly for |^|, I77I, \r\ < 2Lq. Moreover as 
before 

(4.158) the error term in (|4.157|l is 0(iV"") uniformly for |^|, |?7|, \r\ < 2Lo. 

4.3.4. Integrals. Finally, consider the term 
(4.159) 

^Nit V) = 7; f sgn ( r + — - s ) Kn ( s, r + — ) ds, qn as in (|4.138|l. 
Terms of this form arise in the 21 entry of — icl^f'^-'fr + —,r + We want 

qN N,l \ ' ijiv ' qN ! 

to show that lAr(C,??) ^ /q^ ^ Koo{t)dt as N ^ 00. Changing variables s = 
r + [rj + t) / qN we write as 

(4.160) i"jv(e,?7) =Q9,1 + Q9,2 + Q9,3 

where (cf. BTTtIi ^ 



Dr,Nit + ri,T])dt 



1 



Q9,2 = t:I Dr^N{t + r],ri)dt 
1 /""^ 

^9,3 = --;^/ Dr.Nit + r],'n)dt. 
^ Jo 

Again we consider ^, 77, r in the compact set [— 2Lo, 2Lo]- By (|4.150|l . Qg^i converges, 
with all its derivatives, to ^ Kooit + rj — rj) dt — ^ Koc{t) dt as N ^ 00, and 
so we must show that Qg^2 + Q9,3 converges to 0. 
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To analyze Qg, 2,(59, 3 it is convenient to consider the cases |t| < 1 and \t\ > 1 
separately. By the above calculations, we have as ^ 00 

\Dr,N{t + ri,il)dt 



(4.161) 



1 - 

^ A'oo(i + 77 - V) dt + 0(iV-") = 0(7V-"), 



as Koo{t) is even. As before the term 0{N ") is uniform for |?7| < 2Lo- Now 
we must show Qio.i — Qio,2 — > as — > 00, where 

/— 1 /"OO 
Dr.Nit + V,'n)dt and Qio,2 = / Dr^^it + r],rj) dt. 
-co Jl 

Consider (5io,2- We utilize the change of variables 
(4.162) r + {t + T])/qN =CNZ + dN 

i.e. 

t=-T]+ {dN - r)qN + QnCnz 

(4-163) 

z = [r - dMj/cN + (t + r-i)/{qMCN)- 

Fix S > small as in Subsection 14.21 set t± = —i] + (djv — r)qN + qNCN{l ± 6), and 
write Qio,2 = Qua + Qii,2 + Qii,3 where 



Q 



11,3 



Note that t± ~ iV as TV 00. 
We have by (|0?7jl . fTTjl 



;ii,3 



bN-i 



'>N { r+^-^](l)N-i ( r+—]-(l)N-i (r+^-^](f>N (r+ — 



Changing variables t ^ z as in H4.153|l , the first term takes the form 

bN-i4>N{cNZ + div) 4>n~i{cnZo + dN) '^^ ^ 

1+5 qNCNZ + qN(dN - r) - ri 

where zq = z{t = 0) = (r - dN)lcN + vKinCn) = ©(A^^^/^^™)). As b^^i CAr/2, 
using (|3.1|) . H4.16|l we see that the above term is bounded by 

const • {{cN/Nf/^e-"^) ■ [l/cj,'^) < const • e^^^ 

for some c = c((5) > 0. There is a similar estimate for the second term in (5ii,3, 
and we have 

(4.164) |Qii,3l=0(e-^^), c = c(<5)>0, 

uniformly for |r|, \ri\ < 2Lq. 

Changing variables t ^ 2 as in H4.16HII in the first term of Q11.2, and then 
integrating by parts, we obtain 

"1+5 



dt 
T' 



bN~l<pN~licNZo + dN) 



4>n{cnw + c^at) dw 



1 



1-s J 1 + 5 + ^ 

/ 4>n{cnw + dN) dw 



dz 



(z + ^^ 2_)2 

^ Cat a!\rC!\r / 
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which is bounded by 

const • CN ■ c-^"^ ■ {c-^"^N-''^) = 0{N~^'^) 

in view of H4.17|l (of course z ~ 1 in the above integrals, so there is no 

singularity in the integrand). Again there is a similar estimate for the second term 
in Qii,2, and so 

(4.165) |Qii,2| = 0(iV-i/2) 

uniformly for |r|, jryl < 2Lo. 

Let < Of' < min (i, a) and set t^r = iV" . Thus 1 < tjv < L as TV ^ oo. Let 



k2' 



',11 



,1 = Q'ii,i + Qii,i where 



;ii 1 = , Dr,N{t + ri,?])dt and Qn.i = / Dr.N{t + i],r]) dt. 

JtN 

First we consider ^. From 14.143|) we again see that Dr,N{t + tj, rj) is a sum of 
terms of the form S'+ii(- • • )S'+22(- • • ), 'S'+i2(' • • )'5'+2i(' • • ); ' ' ' • Express the term 
S+iiS+22 as a sum 

Q8,i{t + v,v) + QsAt + v,v) + QsA^ + v,v) 

as in H4.144|l . As in the proof of H4.150|) . we have IQsAi^ + ViV)]: IQsait + ViV)] = 
0{N^^) as iV — > CX3 uniformly for 1 < i < and |r|, I77I < 2Lo. Thus the 
contribution of Qs.iit + V,v),Qs.2{t + V,v) to Q'^ is 0{tNN-^) = 0{N-'^+°'') 
uniformly for jryj < 2Lq. From l|4.149|) 

n t-^, \ Q f t + o ( , ^ A sin(7r^G'jv(t + ?7,77)) 

Qs^A^ + ri,V) ^ "3+11 rN H S+22 H 

V qNCN^ V qNCN^ TTt 

where Gn is as in Kim . By Kim 

Qs,3{t + V,V) = S+ii[rN)S+22irN) : ^ i^^^ — I — 

TTt \I\ t>l t 

= S+ii{rN) S+22(rN) cos [n:t{GN{t + r],ri) - 1)] 

COS 7rf 

+ S+nirN) S+22irN) sin [nt{GN{t + 77, 77) - 1)] + 0{N'^). 

But by mSHI), »4.152^ . we find from l|^lT7|l 

GAr(t + 77, ^) = ^^M^l (1 + o(i/7V)) = 1 + o{t/N) + 0{N-^) 

uniformly for l<t<tAr, |77|< 2Lq. Thus 
sin [7ri(GAr(t + 77 , 77) - 1)] 



< \GN{t + 77, 77) - 1| = 0{t/N) + 0(iV-") 



and similarly 

I cos [7ri(GAr(t + r;,77) - 1)] - l| 0(iV^^) + 0(iV^^")- 
Thus as a < 1, 

^ , \ ^ / \ ri ^ \ sinTri / < t 1 ' 

Q8,3(t + ?7, ^) = 5+ii(riv) S+22{rN) ~^^[jp+N^N^^N^ 
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and so 

QsMt + V,v) dt ^S+ii{rN) S+22{rN) / — — dt 



O 



nt 

jy4a' j^2a' j^2a' j^a' 



Assembling the above estimates we find for < a' < min {\,a) as N oo, the 
contribution of the 5+115-1-22 term to Q'^ ^ is given by 



5+ii(rjv) S+22{rN) ^ dt + 0(7V-i+2"') + 0(Ar-«+«') 

uniformly for |r|, \ri\ < 2Lq. Analyzing the other contributions 12 5-1-21, • • • we 
obtain as — > 00 

/*" /■*" sin Trt , I 

Dr,Nit + V,v)dt = J ——dt+OiN-'+^'') + 0{N-^+") 

uniformly for jryj < 2Lq. Again we have used H4.141|l . 
Finally we consider 



Qii.i = I bN~i 



t>N[ r+^—^](f>N-i(r+—]-(j)N~i (r+^-^](j)N (r+ — 



dt 

T' 



Changing variables t — > z as in H4.163|) and then integrating by parts, we obtain for 
the first term in Q'l^ ^ 



bN~i(j)N-i{r + v/Qn) 



1-5 

(pNicNW + diq) dw 



z{tr,) ) QNCNil -5) + {dN - r)qN - T] 

^ ^ {qNCN^dz 



+ / / (t)NicNW + dN)dw . , \ ,9 

Jz(tM) \ Jz(tM) J [INCNZ + (djv - r)qN - vV 

Note that z{tN) ^ N^'^+°' and rccaU (|TT|) . (fH^ to conclude that the first term 
in Qii 1 is bounded by 

const • bN-i ■ c^'/' • (c^'/'a-1) • (^const + QnCn ^) = 0(A""'). 
The same is true for the second term in Q'/i 1, and so together with (|4.166|l we find 
(4.167) 

= ^ Dr,N{t + v,v)dt^ dt + 0(A-i+2°' + A^-"' + 7V-"+"') 

and hence, as a' < 1/2 

Q10.2 = / Dr,N{t + rj,il)dt 

nt 

*" dt + 0(A-i+2"' + A^-"' + A-«+"'). 

TTt 
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The best error estimate is clearly obtained for a' = a* = min (f , I) < ol. Thus 

sin -Ki 

(4.168) y Dr^N{t + r],'n)dt^ J ——dt + 0{N-°') 
and similarly 

(4.169) / DrN{t + V,v)dt^ [ ^^^dt + 0{N-°''). 
J-oo J^tM 

But is even, and together with H4. 161(1 we see that 

(4.170) 

(39,2 + Q9,3 - J ( / ^r,Ar(t + ?/, V) dt - / 7^.,Ar(t + ry, ,7) ) = ©(TV"" ) 



2' ^° 



as — > 00 uniformly for Irj, jr;! < 2Lq. We have proved the following (of. H4.159|l 'l: 
Yn{£., rf) = An{£., 77) + -Bw (??), where as iV oo 
^ ^ (a) yljv(C, 7?) converges with all its 77 derivatives to y Kao{t)dt 

with error 0(A^^") uniformly for ^,ri,r £ [—2Lo,2Lq]; 
(b) Bn{tj) = 0(A^-"*) uniformly for 77, r £ [-2io,2£o]- 
Here a, a* are as in H4.138() . H4.168|l . respectively. 

5. Proofs of Theorems 12.51 12.31 and 12.41 

5.1. Proof of Theorem 12 .51 We know by Lemma ini that D(j)k is a (finite) linear 
combination of 0max(o,fc-n), • ■ • , (f>k+n, n = deg V . Hence 

(5.1) D</.fe= {D(j>k,<l>j)<Pj- 

j— max(0,fc— n) 

Note next that since (j)^ and its derivative are rapidly decaying, 

(5.2) €D(j>k = ^k. ke Z>o. 
Apply e to (|5.1|) and take into account (|5.2(l to find 

/c+n 

(5.3) (j)f. = eD(j)k ^ Y, {D(j)k,(j)j)e(j)j. 

J— inax(0,fc — n) 

Finally, take the inner product of (j5.3|) with any (j)i to find 

(5.4) ^ (Z?(/.fe>,)(e'/'j>i) 

_7— max(0,/c— n) 

which shows that loo — Doo^oo (recall that {D(j)k, 4>j) = for \j — fc| > n). Taking 
the transposes and using the fact that Doc and eoo are both skew symmetric gives 
the second statement in Theorem l2.5r i). 

Now for the proof of Theorem 12. 5r ii). Denote the N x (A^ + n) section of £00 by 
£N,N+n and denote the last n columns of e.N,N+n by 



(5.5) 



B 



12 
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{^*N-nn the same matrix as in ()2.18(l 'l. Denote the {N + n) x N section of 
by Disi-^n,N- Note that the last n rows of this matrix are (0„^jv-n -D2i). By 
Theorem 



IN- 







^N,N+nDN+n,N — ^nDn + ON,N-r 



/„ 

which proves Theorem 12. 5f iil 



B12D21 



5.2. Proof of Theorem I2.3L Recah the definitions of the block matrices A, B, 
and C. Note first that 



BA = 



B12A21 B11A12 \ I — B12-D21 B11D12 



B22A21 B21A12 I \ —B22D21 B21D12 



Recall from Section |2] that R denotes the n x n matrix with all zeros and ones 
on the anti-diagonal. Let denote the transpose of G and be the transpose 
with respect to the anti-diagonal. An application of i? to a matrix from both sides 
interchanges the order of all rows and columns in the matrix. We have for any G 

(5.6) RGR = {G^f = (G^)^ and also RR = /„. 
It is convenient to denote for j, fc = 1, 2 

B,,^7Vi-i/(2")B,fe, D,fc^iV-™2™)i5^.,, i,fe^7V-i+V(2™)A,,. 

Then by Theorem 12. II and 12.21 the matrices with the tilde tend to certain explicit 
limiting matrices, as — * 00, see e.g. (|2.13|) . (|2.3ll . Note that H2.13|l implies that 
as iV ^ 00, even, 

(5.7) = Si2 + o(l), ^^1=^21+0(1) 
and also from the skew symmetry of £00 and D^a 

B12 = --S21, = -D12. 

We prove now the first statement of Theorem 12.31 
(5.8) 

^i?(BA)ni? = RB12D21R = RB^2RRD2iR - {Bi2f(Di^f - {B^^f {D^^f 
= {Bl^ + 0(1)) {bl^ + 0(1)) = {B21 + 0(1)) {b,2 + 0(1)) 

= ^2lL'l2 + 0(1) = S2li:'l2 + 0(1) = {BA)22 + 0(1). 

Next we prove the second statement of Theorem 12.31 First note that 

(S^G)ii {BAG)i2 



BAG 



{BAG)2i {BAG)2 



where 



(5.9) 



{BAC)ii — B12A21 + B12A21B12A21 + B11A12B22A21 
{BAC)i2 = B12A21B11A12 -f B11A12B21A12 
{BAC)2i = B22A21 + B22A21B12A21 + B21A12B22A21 
{BAC)22 = B22A21B11A12 + B21A12B21A12. 
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Let US prove first that {BAC)i2 = 0. Since det A12 ^ (by Theorem 12 . II and since 
det A12 ^ 0) this is the same as proving 

(5.10) B12A21B11 + B11A12B21 = 0. 

We claim that to prove (|5.10() is in turn the same as to prove that B12A21B11 is 
skew symmetric. This is because 

B11A12B21 = BJ^A^^BJ2 = (-812^21-811)^ 

(recall yli2 = A^^). Now since Bn is skew, we note that -612^21-811 is skew if and 
only if 

(/„ + Si2^2l)Bll = {In - Bi2D2l)Bu 

is skew. But the last matrix is skew because it is the lower right n x n corner of 
the skew matrix eNDN^N- (Here we used (|2.18ll and the fact that Bn is the lower 
right n X n corner of e^v-) 

Next we show that {BAC)ii = 0. Since det742i 7^ this is the same as proving 
that 

(5.11) (/„ + Si2^2l)Si2 + -811A12-822 = 0. 

Recall that A12 = -D12, ^21 = ^ -D21 and note that we have to prove 

(5.12) (/„ - Bi2D2i)Bi2 = -B11D12B22 

The matrix (/ — -812-D21) as wc know is the lower right n x n corner of e^DN- In 
view of H2.18|l and (|5.5|) . (/ — -8i2-D2i)i?i2 equals the product of the last n rows of 

[cnDn) with ( I . But this is the same as the product of the last n rows of 

V ^12 / 

Eat with ( ) . We now derive a different expression for the latter product. 

To this end we consider the product of the rows 0, • • • , iV — 1 in Dao with the 
columns N, - ■ ■ , iV + n— 1 in Coo- Since DooEoo = -^00 this product equals ON,n- But 
using the fact that Dqo is banded, and the expression H5.5|l for eN,N+m this relation 
becomes 

-812 / \ — -Di2-822^ 

Multiplying the latter expression from the left by the last n rows of we see that 
the product is —B11D12B22, and so (|5.12|) is established. 

Our next goal is to show {BAC)22 — (-8^)22 + o(l). In view of (|5.8() we must 
prove that {BAC)22 equals 

-R{BA)nR + 0(1) = -RB12A21R + 0(1) = -{RBi2R){RA2iR) + o(l) 
= -{RBi2R){Ai^f + o(l) = ~{RBi2R)Al^ + 0(1) 

= -(i?Si2-R)^12+o(l). 

Now using H5.9|l and inserting tildes, we see that (recall det yli2 7^ 0) we must prove 

-822^21-811 + B21A12B21 + RB12R = 0(1). 

After noting that RB12R = —B21 + o(l), and then taking transposes, we see that 
we must prove 

-812 + -811^12-822 + ^12^21-812 = 0(1). 
But this follows from (|5.11|) (which holds exactly) and which we have already proved 
in the {BAC)ii part. 



D 
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Finally wc show {BAC)2i ~ (-6^)21 + o(l) (which is not really needed in this 
paper; our proof of the fact that liniAr^oo det(/2n + BAC) 7^ in Section |21 only 
uses {BAC)22 {BA)22 + o(l)). By since {BA)2i = B22A21 and again 

det ^21 7^ 0, we must show 

522^21^12 + ^21^12^22 = o(l). 

Since -B22 — -B^ + o(l), this reduces to showing 

B^iA2iBi2 + B2iAi2B^i = 0(1). 

Now we apply the transposition across the anti-diagonal and use the property 
(GH)^ = H^G^ together with (gH) to find 

B11A12B21 + B12A21B11 = 0(1). 

But this follows from H5.10|l (the latter holds exactly) which we proved in the 
{BAC)i2 part. The proof of Theorem 12.31 is now complete. 

5.3. Proof of Theorem 12.41 Assume N is even. Let us consider first the case 
when V{x) is even. Then all the even- numbered OP's and also the functions 4>j 
are even, and all the odd-numbered ones are odd functions on M. As the maps 
/ — > ef, f f reverse parity, it follows that the entries of and 

are zero if i and j have the same parity. In turn we see that the (i,j) entries of 
Wn — ^nDm are zero if i and j have the opposite parity. Let W^^\ ^n'' 
the N/2 X N/2 matrices constructed from the rows and columns 1,3, •• • ,N— 1 
(respectively, 0,2, ■ ■■ ,N - 2) of Wn- Clearly det I^at = det W^^^det But 
more is true: by (|2.18l) and the asymptotics of T^; 2^m-i in H2.16|l . I|2.17|l we 
must have that det M^^^^ = AetT^ + o(l) and deiW'}^^ = detr,„_i + o(l). Now 
Vr^i) = e%^D%^ and W^^^ = (^n^D^n where ej;\ D%\ i = 1,2, are the following 
N/2 X N/2 matrices: 

• ej^\ e^-* are formed from rows 1, 3, • • • , iV — 1 (respectively, 0, 2, • • • , — 2) 
and columns 0, 2, • • • ,N — 2 (respectively, 1, 3, • • • , A^ — 1) of 

• D^^^ , D^^-* are formed from rows 0, 2, • • • , A— 2 (respectively, 1, 3, • • • , A^— 1) 
and columns 1, 3, • • • , A^ — 1 (respectively, 0, 2, • • • ,N~2) of Djy. 

Thus det T^^ = det ej^^ det Zjj^' + o(l) and det T„_i = det ej^^ det fj^' + o(l). But 
as Eat, Dn are skew symmetric, a simple calculation shows that 

Letting A^ — > 00, it follows that detT^ — detTm_i. 

Now for the general case when Vix) is not assumed to be even. Note that the 
A-independent matrices and Tm-i depend only on the degree 2m of V{x) (and 
not on K2m > 0). But this means that detT^ and detTm_i are equal to the 
corresponding determinants for the even weight K2mx'^™'- By the above argument, 
detr^ = detTm_i. This completes the proof of Theorem l2.4l 

6. Differential equations for h{x), ym{0), and proof of Theorem 12.61 

6.1. Plan of the proof. The goal in this Section is to prove that for all m > 2, 
det Tjn-i 7^ 0, where T^-i is defined in (|2.17|l . The key object to understand is I{q) 
in (|2.11(l for q = 3,5, ■ ■ ■ ,2n — 3 = Am — 5. After introducing some convenient nota- 
tion we proceed to the proof. As noted in Remark 12.31 we use different arguments 
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for m "small" and "large," and for a given m, we must estimate I{q) for different 
ranges of q. More precisely, in view of Proposition 16 . II below, what we really need 
is a bound on |1 + i{q)\, where i{q) = ml{q) — \ (see (|2.20|l and (|6.3() below). For 
some values of q and to, we estimate |/('z)|, and hence |l + /(g)| < 1 + |/(q)|, but in 
other ranges we must consider the absolute value of the combination 1 + I{q). Of 
course if we could show that 

(6.1) \i{q)\ < 1, 

for any to > 2 and all g = 3, 5, • • • , Am — 5, then Theorem 12 . 61 would follow immedi- 
ately from Proposition 16. II unfortunately we are unable to prove (|6.1I) for all q,m 
as above, and in fact it may not even be true for all q, m in the range. 

In Subsection 1^21 we prove a crude a priori bound on the absolute value of I{q) 
(see (|6.15|) below) that is valid for all to > 2 and all g = 3, 5, • • • , 4to — 5. The proof 
of this bound follows from the hypergeometric ODE for h{x) mentioned above (cf. 
below). 

In Subsection l6.4l we use a more refined argument to show that the absolute value 
of I{q) is bounded by 1 for all to > 2 and for all q in a region 0{y/rn) < q < 0{m), 
see H6.21|l below. To this end we use the properties of ym{d) (see H2.21|l ^ that 
follow from the fact that it satisfies a certain Riccati equation (see Proposition 
16.4(1 : as noted in Remark 12 . 31 this Riccati equation in turn is a consequence of the 
hypergeometric ODE for h{x) mentioned above. 

Finally, in Subsection 16.51 we use a certain integral representation for h(x) to 
find, for large enough m, an accurate approximation for I{q) which yields a bound 
(see (|6.57(l below) on |1 + I{q)\ for q in the region 3 < q < 0{y/m). 

The above estimates are sufficient to prove det Tm^i ^ for all to > 2. The crude 
bound in Subsection 16. 31 is sufficient to give a simple proof of the result for to < 51 
(see Subsection l6.3|) . The bound in Subsection l6.4l for 0(-/to) < q < 0{m), together 
with the crude bound in Subsection 16.31 applied to q in the range 3 < 9 < 0{y/m), 
gives a proof of the result for to < 99 (see Subsection 16.411 . Finally the bound 
in Subsection 16.51 for large to and 3 < q < 0{^/m), together with the bound in 
Subsection 16.41 for 0{^/m) < q < 0{m), proves detTm-i 7^ for to > 38 (see 
Subsection 16. 5|l . In this way the estimates in Subsections 16 . 31 and cover the 
entire range to > 2, with overlap and some redundancy. 

6.2. Notation. Recall the definition of h{x) in (|2.8|) . We will repeatedly use the 
change of variable x — sin 6 in (|2.11() . We find 

(6.2) I{q) = + A sin ^ / ^ cos{q0) y„ {9) dO 

2m TOTT 2 

where ym is defined in 12.21|l . and as indicated above we introduce the notation 

(6.3) /(g) = TO/(<7)-i = -sin^ / ^ cos{q9) y„,{e) dO, g = 3, 5, • • •4to-5. 
In H6.2|) we have used the following, simple identities for q in the above range, 

(6.4) cos{q0) cos 0d0 = O, - sin ^ /"^' d9 = 2. 

^ ' J-./2 ' TT 2 COS0 
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Set 



n2 



2(m!) 

Im = TZ 77 

myzmy. 

and let i^m-i = iXjkYj'kli denote the very last matrix (of binomial coefficients) in 
(|TT7jl . Set 



( ' ' 

1 1 
V 1 1 



1\ / i{n) i{n + 2) 
1 ^ /(n-2) i{n) 

V V /(3) /(5) 



/(2n- 3)\ 
I{2n - 5) 

i{n) ) 



Then (|TT7|) becomes 
(6.5) 



We treat the second term as a perturbation and it turns out to be important to 



choose the norm 

pm — 1 



on appropriately. We take || • || to be the maximum 

norm on M'"^", \\x = (xi,--- ,a;m-i)|| = niaxi<j<m,-i ja^j |- We show eventually 
that acting as an operator on (K.™^^, || • ||), the second term in Ht).5|) has norm < 1 
for all TO. This of course implies that detT,„_i ^ for all m > 2. 

Proposition 6.1. Assume that for some C > 0, maxj.fe^i.... .m-i ^ C. Then 



(6.6) 



|| S ^77 7Z 77 ' ^ 



.2 to(2to)! 
Proof. By the definition of the norm || • || on R™~^ 



l-'^^m-i^m-ill < max > \{X„i-iY„i^i] 

2— I,--- ,m — 1 ^ — ^ 



m— 1 m— 1 

< . max ^ <C ^ 

z— I.--- .m— 1 ^ — ' — ^ 

i,k=l j,k=l 

m-2m-2~i ,^ s m-1 i\ 



1=0 ;=o 



Substituting X]/!lo 



m-l (2m-\\ 



= 2^ 



and 



2 V 1 



which follows from the identity 



(6.7) 



i=0 



9-1 



9 < 
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with q = m — 1, p — 2m — 1, we arrive at H6.6(l . It remains to prove H6.7|l . Note 
that 

1=0 ^ ^ 1=1 ^ >' 

9-1 / \ 9-1 



9 / \ 9-i / \ 9-i / \ 

1 = 1 ^ ^ 1=0 ^ ^ 1=0 ^ 



from which H6.7() follows. □ 

6.3. ODE for h{x). The final estimate (I6.16|) in this Subsection is completely 
superceded by the final estimate (|6.22f) in the next Subsection, which gives the 
desired result for m up till 99. However both H6.16|l and H6.22|) utilize the estimate 
(I6.15|l that is proved below in this Subsection. 

Proposition 6.2. For any m G N the (even) function h{x) satisfies the ODE 

(6.8) x{x'^ - l)h' + {{2m - 1) - 2{m - l)x^)h = im, 0<x<l 
Also 

(6.9) h{Q)^-^, Ml) = 4m 

2m — 1 

which together with 16.8|l imph 



ICS 



Am /"^ /a;\2"i dt 



< x < 1. 



1 I nmn n li^ i Q ^ 

2m- 1 



Finally, the function h{x) is strictly increasing on [0, 1], and hence it is > r^^_^ > 
on [-1,1]. 

Proof. Assuming (|6.8|) and either of the conditions (|6.9|l the integral representation 
(|6.10|l can be found using the integrating factor method. The explicit formula (|2.8() 
readily implies H6.9|l after elementary algebraic computations, and also that h{x) is 
strictly increasing on [0, 1]. 

To prove note that from lfTH|) 

_ 4to (-TO+l)fc 

~ 2m- 1 (-TO + 3/2)fc 
where {a)k = a{a+l) ■ ■ ■ {a + k— 1) denotes the Pochhammer symbol. As (l)fe — kl, 



4m {l),{-ni+l),{xy 

n — 

Am 



h{x) = V 

(6.11) ^^-^h (— + 3/2). fc! 



-2Fi(l,-m+ l,-TO + 3/2;.T^) 



2m- 1" 

(the series in 2^1 terminates at fc = m — 1). Thus w{z) = h{^) satisfies the 
hypergeometric equation 

z(l - z)w" + (c - (a + 6 + l)z)w' - ahw = 

with a = 1, 6 — — m + 1, c = — m + 3/2 which becomes 

(6.12) z(l - z)w" + ((-m + 3/2) + (to - 3)z)w' + (to - l)w = 0. 
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= 



In turn this equation can be rewritten as 

(6.13) ^ (z(z - l)w' + ((m - 1/2) - (to - l)z)n?j 
and so 

(6.14) z{z - l)w' + ((to - 1/2) - (to - l)^)^ = const. 

Letting z ^ 1 and using (|6.9(l we conclude that const = 2to. But H6.14|l is just 
(|6.8(l for w{z) = h{^/z). Of course 1)6. 8|l can also be proved by substituting H2.8|l 
directly into (|^ . □ 

Remark 6.1. Differentiating H4.i;-{l|l and using H6.8|l . we immediately obtain the 
proof of Lemma [4.81 

Now from 16.8|l it follows that 



1 2(TO-l)-^^'(^) 



1 — \ h{x) I h(x) 
which together with (|6.3|) and (|2.21|) gives for q = 3, 5, • • • 

~ TO 2 qn cos(q arcsina;) / 4m \ 

/(g) — sin — / , , ^ — I] dx 

4TO7r 2 i„i l-a;2 \h{x) J 

1 , qn . N /„/ X xh'(x)\ , 

= — sm — / cos((7 arcsmx 2(to — 1) ——]ax. 

27r 2 i_i ^ ^ \ ^ ' h{x) J 

The first integral vanishes for odd <? > 3 by (|6.4|l . and estimating the second by the 
absolute value of the integrand (using ft,' > on [0, 1] and the fact that h is even) 
we obtain the a priori estimate 

1 xh'(x) , 1 /"W, h(x)\' 1 f\ 4to , 

71" Jo Hx) TTJo V h{l)/ TT Jq h(x) 

by (|6.9|) . Now we use the estimate 

— - < (2to- 1) -2(to- l)a;^ < a; < 1, 
h[x) 

which follows from (|6.8|) since h, h' > 0, to find 
(6.15) \iiq)\<L{m) 
where 

_ l0g(2TO - 1) - 2 + {t,n + 1) log(l + t-,1) - {t„, - 1) log(l - t-.l) 



i(m) 



TT 



and tm = y 2m-2 - Using Maple, for example, it is straightforward to check that 
for TO = 2, 3, • • • ,51 

(6.16) (1 + L(to)) f i - • 22^-2) < 1 

^ ' ^ ^ \2 to(2to)! / 

and hence by Proposition 16.11 T,„_i is invertiblc and therefore detT,„_i ^ for 
TO < 51. 

Remark 6.2. By Stirling's formula, i^'J,'; ■22'"-2 ^ ^nd so the LHS of 

grows logarithmically as to — )■ oo. 
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6.4. Riccati equation for ym{0). A consequence of the calculations in this Sub- 
section is that 

l!7m-'^m-iym-l|l(R'"-M|-||) — (R^-MHI) < 1 

for m up to (the larger value) 99. In place of Proposition 16.11 we will use the 
following result. The proof is immediate. 

Proposition 6.3. Assume that for some matrix X with nonnegative elements we 
have for all j, k — I, ■ ■ ■ ,m — 1 that \Xj^k\ < ^j.k- Then 

m— 1 

(6.17) \hrnX„i-lY„i-l\\ < Jm Hiax XikYkj . 

.m— 1 ^ — ^ 

j,k=l 

The main goal in this Subsection is to prove the estimate (|6.21() below on I{q). 
This estimate is useful for q in the range 0(-\/to) < q < 0{m). The matrix X in 
Proposition 16.31 is then constructed by "interpolating" between this estimate and 
the estimate H6.15|l above, as indicated in H6.22(l below. Estimate (16.21(1 will be 
derived from the Riccati equation ^121 for yrn{0) = ^(^^Te) " - ^) 
given in Proposition 16.41 below. Note from (|2.8|1 that ym{9) is even. 

Remark 6.3. The specific form of ym{0) has the following consequences for /(g) = 
|-sin^ ^"^^ '^^'^ show that as m — *■ oo, ^^^T~ ~ 

(formally at least, divide (|6.8|l by m and let m oo). Hence for a fixed q, as 
m oo, we expect I{q) — i (recall (|6.4|) 'l. On the other hand for odd q ^ 21 + 1, 

cos{qe)ym{0) = cos(2/0) (cos0 2/™(0)) - sin(2/^) {sine yrn{Oj) 

and so I{q) is a sum of Fourier coefficients of the real analytic 7r-periodic functions 
(cos 0) 2/m(0), (sm9) yjn{0), and so for a fixed to, I{q) decays exponentially as q — > oo 
by the Paley- Wiener theorem. As discussed in Remark 12.31 above, this indicates 
that there should be a transition in the (g, TO)-plane between these two kinds of 
behavior. And indeed, as we will see below, the transition region is in the range 
q ^ yfm. 

Proposition 6.4. For any to > 1 the funetion ym{0) defined in H2.21|l satisfies a 
Riccati differential equation 12.2211 

4 / 2TO+1 \f I \ 

Vm = -T-^ Vrn + -. COS 6* y„ + . 

sm \ 4 / \ 2 cos J 

Also yrniO) ~ —1/2 and j/„j(7r/2) = 0, see Fig.^ 

Proof. This is a straightforward computation using the ODE (|6.8|l for h{x). The 
interesting fact here is that the RHS in the Riccati equation has real roots which, 
moreover, have simple expressions. The proof that ?/m(0) = ^ 1/2 follows di- 
rectly from H2.9|l and the fact that ym{T^/'2) = follows from H2.9|l together with 
L'Hopital's rule. □ 

We will show next that y„i is a unimodal nonpositive function on [0, 7r/2], and 
that its minimum satisfies a certain bound (see Fig.|21). 

Proposition 6.5. For any m > 2 the function ym{0) on the interval [0, 7r/2] strictly 
decreases from ym{0) ~ —1/2 to the value ?/m,min which is achieved at a unique point 
Omin G (0, 7r/2), and then strictly increases to at the point — 7r/2. The value 



UNIVERSALITY FOR ORTHOGONAL AND SYMPLECTIC ENSEMBLES 69 



0.2 0.4 0.6 0.8 1 1.2 1.4 




Figure 1. Graphs of ym{d), ^ cos6> for ni — 10 and 
-1/(2 cos 6*) (the latter for < < f - 0.1) 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 




Figure 2. Graphs of y„^(B) for m = 10, 100, 1000 and - j^i 
(the latter for < 6* < f - 0.04) 



Vm.-miw lies above the y-coordinate of the point of intersection of the curves — ^ ^ 
and — ^"^+^ cos 9, see Fig.^ Finally, ym{S) < and 

(6.18) ym,min>-\^f^^Tlj2. 

Proof. It is helpful to refer to Fig. ^ in the following argument. Define u{9) = 
and lynlo) = -^2|hicos6l. Note that by (ESI) and itT^ it follows that 
ym{9) is a smooth even function on [—tt/2, 7r/2]. This implies y'mi^) = 0. Also it is 
clear that for aU m > 2, y,„(0) = u(0) > ;,„(0). Also m'(0) = and it"(0) = -1/2. 
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Assume for a moment that 

„ , , 1 2m — 1 
(6.19) y-^'^--22^3- 

Then ?/''j(0) < —1/2 for m > 2. Hence for small 9 > 0, ym lies below u and above 
Im- It follows in particular as ?/m(7r/2) = and u{0) -co as 7r/2, that i/m 
must intersect u at some point, say, in (0, 7r/2). 

It the argument that follows the signature table for the vector field 

sm a 

plays a crucial role. We have 

• y{ym, 6*) > if ym hes above Imid) and u{9) or below /„i(^?) and u{9) 

• y{ym, 6*) < if ?/,„ lies between and w(6'). 

Let 6 E (0,7r/2) be the (unique) intersection point for lm{9) and u{9), lm{9) = 
u{9) = y. From the signature table we see that V{y, 9) < for all 9^9. Moreover 
if y^i9) = y then y'^{9) = V{ym{9),9) = 0. But u'{9) < and C(0) > 0, and 
so ym{9) lies between lm{9) and u{9) for 9 > 9, 9 close to 9. From the signature 
table it then follows that ym{9) decreases as 9 crosses 9. In other words, we see 
that if ym{9) crosses the level y, then it must decrease. As 2/m(0) = —1/2 > y, 
and as ymij^/'i) = 0, it follows, therefore, that ym{Q) cannot cross y. In particular, 
Q < 9q < 9 . Now a similar argument shows that for 9 > 9o, 9 close to 6*0, ym{9) 
lies in the region above lm{9) and u{0) (and increases across ^o)- Suppose that 
ym{9) exits this region for the first time for some 9i, 9q < 9i < tt/2. But then, 
again by the signature table, 2/^(0) = V{ym{9),9) > ior 9o < 9 < 9i, and 
hence ym {9i) > ym{9o)- In particular, as u'{9) < 0, it follows that ym{9) does 
not exit the region through the curve u{9). Thus we must have ym{9i) = lm{9i) 
and 9 < 9i < 7t/2. But arguing as before we see that ym{9) must decrease as it 
crosses 9i: as ym(7r/2) — there must be some point 92, 9i < 92 < Tr/2 for which 
ym{92) = y7n{9i) and y'jn{d2) > 0. But F(y,„(6'2), 6*2) < 0, by the signature table, 
which is a contradiction. The above arguments show that y„i(0) crosses u{9) at a 
unique point, 9nwi — 9q, < 9o < 9: for < < 9o, ym{9) lies between lm{9) 
and u{9) and for 9o < 9 < tt/2, ym{9) lies above lm{9) and u{9). Thus ym{9) is 
unimodal, decreasing for < 6* < 0o and increasing for 9o < 9 < tt/2. Clearly 
2/min — ?/m(^min) hss abovc y = u{9) ~ lm{9) ~ ~\\J m + 1/2. Note finally that for 
< < ^niin we certainly have ym{9) < u{9) < 0. On the other hand as y'„i{9) > 
for (?niin < 9 < tt/2, and as ym{TT/2) = 0, it follows that we must have ym{9) < 
for all 6'min <9 < tt/2. Thus y„(6') < for all < 61 < 7r/2. 

It remains to prove (|5.19|) . Divide the ODE (|2.22|) by sin 9 to obtain 



smt/ V 4 y sin^l 



2cose 



Now as y„,(0) = -1/2, 

e^o sin^ 2 sin 6* cos 6* 2 6^0 sin0 4' 

Letting 61 -> on both sides of (jOUjl we obtain hm^^o = -^f^- But 

7/^(0) = linie^o ^^^tI^j and hence (|6.19|) follows. This completes the proof of the 
Proposition. □ 
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1/(9)1 



< 



r/2 

97r \ Jo 



4 









r/2 



sin(g^?)y:„(0)d0 

/2 



where we have used y„i(7r/2) = and also the unimodaUty of fProDOsition l6.5|l . 
Hence in view of (16.18(1 for any m > 2 and any odd q > 3 we have 



(6.21) 



In particular, \i{q)\ < 1 for odd q such that + l/2j + 1 < g < 4m — 5. 

Now we can combine the estimates (|6.15ll and ()6.21|l at our discretion. We define 
the estimating Tocplitz matrix Xm-i = {^j,k)J^k=i Proposition 16.31 as follows. 
The diagonals are numbered by odd q = 3,5, • ■ ■ , 4m — 5 from the lower left corner. 
Thus if j, fc = 1, 2, • • • , m — 1 are the row and column indices starting from the 
usual upper left corner, then q — 2{k — j) + (2m — 1). Therefore, in view of I(fi.l5|l . 
(|OT|l . set 



(6.22) 



X 



3,k 



(2(fc-j) + (2m-l))7r 



y/m+1/2 



Also let Yrn-i = {Yj,k)J^^2i where Y,- & = C'^-/) for 1 < j < A: < m - 1 and 
otherwise, be the matrix of binomial coefficients as before. Then one can check 
by computer (only sums and products are involved) that the RHS in (|6.17|l is < 1 
for m = 2,3, •• • ,99. Thus by Proposition 16.31 r,„_i = /,„_! — ^mXm-iYm-i is 
invertible for m < 99, and hence detT^-i ^ for m < 99. 



6.5. Large m asymptotics for ym{d) via an integral representation. This 
is the last and most technical computation. The main idea here is to find for 
large enough m an approximation to the solution ym{G) of the Riccati equation 
in J^^'^ cos{q9)ym{0)d0 such that the integral is well approximated in the range 
5 = 3, 5, • • • , 0(-\/m). We show that the estimate (|6.57|l below that follows from 
this approximation, together with (|6.21|) for [|--\/m + l/2j +1 < q < 4to — 5, gives 
the desired result for all m > 38, so that there is a large overlap with the region 
m < 99 where the result of the preceding Subsection is valid. 
For p > set 



(6.23) 



G{p) 



/o\l+p«)-™-l/2. 



dv 
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The first result we need is the following integral representation for ym{0)- Using 
X = sin 9, 

... , 1 1 l^x^- 

ym{0) 



Vl - a;2 \h{x) Am 



(6.24) = -^lEZ _ -^(m - i)^^'(^)'"((r-l)7fe 



Vl-^a? rn — 1 



where 



2 

cos^e* 



- x2 G{p) 



(6.25) p=a;"^-l = — x^sin^. 

sin 

We will explicitly indicate the dependence of h{x) on m as hm{x) when needed. To 
prove (|6.24|l we first note that for m > 2 

2m 

(6.26) h„,{x)^- -(2 + a;2/i„_i(x)), < x < 1, 

2m — 1 

which follows from the explicit formula H2.8() and in turn leads (after adding and 
subtracting 2^-1 ('^ ~ l)^m) to 

(1 + (2m - 2)(1 - x^))hm = 4m + 2x^{mhm-i - (m - l)/im) 
from which, dividing by Amhm, we obtain 

1 1 1 — .T^ 1 — x^ x"^ mhra-i — {m — \)hj,^ 
h„i Am 2 2m 2m hm 

Together with (|6.10() , this implies the second equality in H6.24|l . The third equality 
in (|6.24|l follows after changing the variable t — x^l + pv where p is as in H6.25|l . 
Now from 16.24() for g = 3, 5, • • • , and H6.4|l . 

m — 1 OTT f^/^ 

I{<1) = 2sin-^ / { cos{q9) cose) G{p{9))d9. 

Make now a change of variable cf) = j — 9 giving 

(6.27) p = tan2(/), dp = cos =—=!=, sin(/)=^^= 

cos-' V 1 + P V 1 + P 

to find 



TT Jo ^"^(l + p)3/2- 

We split the integral into two, 

(6.28) i{q)^-!!l—( + ) sin(g0(p)) G(p) ^ 



(l + p)3/2' 

Assume that m > 3 (we only need the results below for m > 38). 

Lemma 6.6. We have 

m-l sm{q(t>{p)) G{p) m-1 1 p sin(#(p)) 
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where the error term Ei{m,q) satisfies 

Denom[m) 

uniformly for all q ^ 3,5, ■ ■ ■ , where 



A{m) 



7 {m-l)y/m 



ttV2 (m- 5)(m- |) 
B(m) = ^ -(logm) TTTT + - 



(6.29) 



C(m) 



V-^ ^ 2 
1 



V \/2 

Denom[rn) = f-l 1 — e . 

m- 2 V / 

Proof. We consider p > to^^/^. We need the following elementary inequality: 

(6.30) e'^y/^^+y^ <{l + yY <e°'y, a,y>0. 
Set 

(6.31) fc = m + l/2. 

Splitting the interval [0, 1] into two and integrating by parts we find 



(l+pv) 



dv _ 1 V2 1 



(fc-i)p (fc- i)p (i + p/2)*=-i 

1/2 



and 
(6.33) 



1/2 Vl - 1' 



, ,1, vdv 1 111 

(l + pv) 



(fc-l)(fc-2)p2 (fc-l)p V2 (l + p/2)'=-i 

1 3 1 

" (fc-l)(fc-2)p2 7f (l + p/2)fe-2 



1 



1/2 



-k 



1 1/2 \/l-V 

of the ratio in G{p) for ' 
1 (fc- 2)p- (RHS in iniSSl) - (R-HS in in!S3l) 



This gives the leading behavior of the ratio in G{p) for m (and k) large and motivates 
the definition 

A{p) ^ G{p) 



ik-2)p (fc_2)p/o(l + pi;)-fe-7fS 



l-u 
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which yields 



1 k-2 1 



(6.34) 



(fc-l)p fc - 1 (1 + p/2)fc-i 

1 3 1 
~ (fc - l)p 71 (1 + p/2)'=-2 



+ (fc - 2)p / (1 + p^;) 



-fc 



vdv 



'1/2 \/l - 

1 V2 1 



(fc - l)p (fc - l)p (l + p/2)'^-i 



(fc 



dv 

/2 



From this we derive an upper and a lower estimate on A(p). Discarding the negative 
terms, we obtain 

((fc-2)p^\l + pt;)-'=-^) A(p) 

(6.35) (fc-l)pio ^ Vx/T^; 

+ {k-2)p[ {l+pv)-^ 



+ 



1/2 \/l - 1^ 

y/2 1 



(fc-l)p (l+p/2) 



fe-1 ■ 



Note that the second derivative above is indeed nonncgativc and that its maximum 
over [0, 1/2] equals 7 /\/2. Next changing variables \/l — v = t and integrating by 
parts we obtain 

(1 + pvY^ , =2 / :r—rdt<2 —r:dt 

1/2^ ' 7o (i + p(i-t2))fc - 7o (i + p(i_f))fe 

(l+p(l-l/V2))~^+i l 2 /■^/'^ -2t 

(fc-l)p 2 (fc- 1)^70 (l + p(l-t))'=-i 

1 2 ^1 /"^/^ 



(fc-l)p(l + p(l 



2 /I /"^'^^ \ 

^ -+/ 2tdt\. 

l/V2))*=-i V2 io ) 
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Thus from ^^B, 
(6.36) 

f (m - 3/2)p / (1 + H"™~'^' -^^) < ^7 77:^77^^ fTTTTT 

' ^'^Jo yr^/ ^'^^ - \/2(to- l/2)(m-5/2)p2 

2 ^/2 1 



+ 



(1 + (1 - l/V2)p)™-i/2 (to - l/2)p (1 + p/2)"-i/2 ■ 
Next discarding the positive terms in (|6.34|) . we find 



-((fc-2V^\l + H-'^^) A(p) 

fc- 2 1 1 13 1 

(6.37) < 



A: - 1 V2 (1 + p/2)^-i (fc - l)p V2 (1 + p/2)'=-2 



{k~l)pj„ V\/l - 1'/ A/2 \/T^ 



Note that the first derivative above is indeed nonegative and that its maximum 
over [0, 1] equals V2. Note finaUy that J^^^{1 + pv)-"^ < ^f+jj^ Jy^ ^ = 



Thus (|n23, |OT|) give 



(l+p/2)'= 



-((m-3/2)p Al + p«) — 1/2-^) A(p) 



(6.38) < 



1 TO- 3/2 1 3 1 1 

V2 m - 1/2 (1 + p/2)'"-i/2 + V2 (m - l/2)p (1 + p/2)"-3/2 

\/2 V2 
(to- l/2)(m-3/2)p2 + (1 +p/2)™+i/2- 

Now using H6.3()|l . we see that the factor in front of A(p) in H6.36|l satisfies 
(6.39) 

{l + pv)-''-^> f\l+pv)-'^dv = —^(l- ^ 



- {k-l)p{ J - (fc-l)p' 

> 1-e ^ 

- (to-1/2)pV 

by (|6.31() since —j^ decreases and to^^^^ < p 

We now consider the various terms in H6.36|l . H6.38|l to obtain an estimate on 
|A(p)|, p > to1/2. Note that: 

• the 1st term in H6.36|l dominates the 3rd term in (|6.38|l : 

• the 2nd term in (|6.38(l dominates the 3rd term in l|6.36|l : 

• both the 2nd term in H6.36|l and the sum of the 1st and 4th terms in 

are dominated by ^.^^.J/j^^^^^.^,, - 



76 



DEIFT AND GIOEV 



Thus taking into account H6.39|l and recalling the notation (|6.29(l we conclude 
(6.40) 

,A/ M 1/7 1 

|A(p)| < 



Denom{m) \y/2 {m — l/2)(m — 5/2)p^ 

3 1 13 1 



(m - l/2)p (1 + p/2)™-3/2 V2 (1 + (1 - l/V2)p)™-i/2^ 

for all TO~^/2 < p < oo. 

Now substituting G{p) = (^_3/2)p + -^Ip) ii^to the integral in H(j.28|l . we 
obtain 

m-1 /-"^ sin(g0(p)) G(p) m-1 1 sin(g0(p)) ^ , ^ ^ 



1/2 (l+p)3/2 ^ m - 3/2 TT 7„,-i/2 p(l + p)3/2 



TT 

where 
(6.41) 

1 m - 1 / 7 1 f°° dp 



Denom{m) tt ^2 (™ - 1/2)(to - 5/2) 7„-i/2 _|_ p)3/2 
3 1 f°° dp 



V2 (m - 1/2) 7„,-v2 p{l + p/2)™-3/2(l + p)3/2 
3 f°° dp 



V2 (1 + (1 - 1/V2)p)™-1/2(1 + p)3/2 , 

After elementary manipulations (splitting the integral as below) the first term gives 
A{m) in (|6.29|) . To estimate the second term, note that 
r-l roo . 



/2 7i yp(l + p/2)"-3/2(l + p)3/2 

1 dp 1 /'°° d/9 



< 



(1 + 1/(2V^))'"-3/2 p (3/2)^-3/2 7^ ^5/2 

which gives B{m) in (|6.29|) . Finally in the third term in H6.41f) we use 

dp ^ f°° dp 



-1/2 (1 + (1 - 1/V2)/9)'»-1/2(1 + p)3/2 - (1 + (1 - 1/V2)p)"+1 

and arrive at C{m) in 16.29|l . This completes the proof of Lemma IH?^ □ 
^-1/2 ____ 

In the integral in H6.28|l we approximate the function G{p) by the function 

(6.42) Ga{p) = A 

Jo ^ 

and use the following result. 

Lemma 6.7. For m > 38, we /iat;e 
(6.43) 

m - 1 r sin(#(p)) Gip) m-1 r'"'' sin(g0(p)) ^ . . , ^ . . 

dp^ / — — — ——G^[p)dp+E2{m,q) 



(l + p)3/2 ^ ^ 7o (l + p)3/2 

where the error term E2{m, q) satisfies, uniformly for all q ~ 3,5 

\E2{m,q)\ < -^-6, 
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where the number 6 provides the following (crude) bound (proved in the Appendix) 

1 ~-svR{m) c 
)_ 

7n>38 s>0 „-sv dv 

Jo ^ v^T^ 



(6.44) max max s^-^^ — j < 6, 



where 



(6.45) R{m) ^ YT^- 

\/ni 

As R{m) increases with m, the maximum on the LHS of ()6.44|l is achieved at 
m = 38. 

Proof. We consider < p < m^^/^. Recall the notation (|6.31|) and denote 

^k[s) = S -J . 



e' 



-sv _ 



Assume for a moment the estimates 

(6.46) o<G(p)-G,(p)< 0<p<m-i/2. 

Then (|6.43|) holds where E2{m,q) satisfies 

-1/2 

M m — 1 1 f™ . ,-, . dp 



n kJo " (l + p)3/2 

-1/2 



/ \ 1 /" do 

< ( maxAfc(s))- / ,^ , ,,,, <6 

V s>0 ' ]q (l + p)3/2 



'^p ^ ^ 1 1 



uniformly for all m > 38 and all g = 3, 5, • • • . This is the desired result. 

It remains to prove (|6.46|) . First consider the upper bound. In view of H6.23|l . 

G(p) < ^° = G.(p) + ^° ^ ^ ^ 



„~-kpv dv ""^ ' „-kov dv 

Jo v/T^ Jo 



Next, for aU < w < 1 

g-fcpi,/(l+p) _ g-fcp« ^ (^-kpv{l~-t+t/{l+p))\ ^ ^P^" /"^ g-fep«(l-t+t/(l+p)) 

where the exponent is estimated by its maximal value at t = 1 and we use the fact 
that —ji^ < —R{m) for < p < by the definition H6.45|l . This proves the 
upper bound in H6.46|l . 
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For the lower bound in H6.46|l we note 



{G{p) - Gaip)) • ( / + -^=^) ■ ( 



-kpv 



dv 



[(1 + pi;)-'=^e~'=''"' - (1 + pw)-''e-^P''j {v - w) 



dv 



dw 



"'0 
1 rl 



V ^/l — w 



2 Jo Jo 



{v - w) {1 + pv)-'' e-'^P" 



1 + pv 
1 + pw 



dv 



dw 



\/\ — V vT— w 



Consider J{x) = x^ fc > 1. Note that for a, 6 > 0, /(a) - J{h) ^ f{e{a, b)){a - h) 
for some 9{a,b) between a and b. As 9{a,b) > 0, we must have f'{d{a,b)) > 0. 
Thus 



(6.47) 



{G{p) - Gaip)) • ( / '(1 + Pv)-' ■ ( 



-kpv 



dv 



Jo 



\/l - Wo ^/l-v ' 

dv dw 



\v~w) (eP^^--^ - F,Av, w) 

\ 1 + pw J 



w 



where Fk^p{v,w) > everywhere. Consider now fwiv) = e^'^^ — jj^- Then 
U{w) = and /^(w) = p{eP'-'"-'^^ - j-^) > for w > w > (note p > 0). Thus 
Uiv) >OfoTv>w. Also Uiv) = e"'"— )i±£^(i±^ - eP(— )) which is < by 
the same argument ioi < v < w. We conclude that the indegrand in l().47|l is > 
for all < w < 1 which proves the lower bound in l|6.46|l . The proof of Lemma 
16.71 is now complete. □ 



The integral that appears on the RHS in (|6.43|) . may be written in the form 



m-1 sin(#(p)) 

(6.48) {{m + ll2)pGM - l] dp 

m - 1 1 r'"^ sin(#(p)) , 
' dp 



m+^TTJo p(l + p)3/2 
The second integral on the right in H6.48|l equals 

m - 11 r~"'' sin(g<^(p)) m - 1 1 f''' sin(g^(p)) 

■dp = 3-/ — — — —dp + E3{m;q) 



m 



^nJo p(l + p)3/2 m- I Trio p(l + p)^/ 



UNIVERSALITY FOR ORTHOGONAL AND SYMPLECTIC ENSEMBLES 



79 



where 



„_2 m-1 p"'' |sin(#(p))| ^ 

\Es{m;q)\^-j-^^^^ 3^^^ -7^—^ dp 

-1/2 



'tt (m+ i)(TO- |) 7o p(l + p)3/2 

2 g m-1 \(b{p)\ 

p(l+p)3/2 

-1/2 



< 



< 



TT m + 


1 
2 


m — 


3 
2 


2 9 




m — 


1 


TT m + 


1 
2 


m — 


3 
2 


2 g 




m — 


1 



dp 



(6-49) < Z-^^ / -77^^ 



-1/2 



p(l+p)3/2 

dp 



4 g m — 1 1 



TT I 1 m — ^ m 



3 3/4 



and we have used 

(6.50) < tan(0(p)) = 

by (|6.27|) . Denote the first integral on the right in (|6.48B by £^4(m; q) and set 

Jo ^ 

Then (m + l/2)pGa(/9) - 1 = F{{m + l/2)p) and 



I £'4 (to; g) I 



"^-^^ ^""^ ^^^('^'^(^)).F((TO + l/2)p)dp 



m + ^n Jo p{l + p)3/2 



<-g / |F((to+1/2)p)|^ 
Jq ^/P 

1 q '•+°° 



< --^ / 2\Fis^)\ds 



71" /™ I 1 







where we have again used (16.50(1 . Now we denote 



1 ^s^v V dv 



lo e 



g(.).2n.-)^2(,s- - 
Jo ^ 



and refer to the (crude) bound proved in the Appendix 

(6.51) / \H{s)\ds< 2.8 

Jo 
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to conclude that 

(6.52) |^4(TO;g)| < - , ^ ■ 2.8. 



Combining Lemma 16.61 and Ifi . 71 we find 

m - 1 1 f°^' sm{q(f){p)) 



4 



Returning to the variable 6 = ^ — (t>{p) and using (|6.27|l we see that for any 
1 r-(#(p))^^^2 r/.3in(^^^^^_ 



7i"7o p{1 + pY/^ it Jo sin 
Hence (|6.53f) becomes 

5 

(6.54) /(g) = -l + Ei?,(m;g) 
where 

11 ,,,,11 

E5[m-q) = -- 3-, \E5{m;q)\<- 3. 

Zm— J 2 m — ^ 

Recall that we need a bound on 1 + /(g). Assembling the above estimates for 
m > 38, g = 3, 5, • • • , Am — 5, 

5 5 

1 1 + i{q) I = I E ("^; 9) < E 1^ I 

(6.55) ^ A(m) + i3(m) + C(m) ^ 16^1 1 

~ Denom{m) tt 2 m — | 

g /2.8 4 m - 1 1 



+ 1/2 V 71" TT m - I 7713/4 ^ 

We can now complete the proof of Theorem 12.61 for the case of large 777. We 
wish to apply Proposition 16.11 We estimate the elements of the matrix Xm^i as 
follows. Enumerate the diagonals of Xm-i by g = 3, 5, • • • , 477i — 5 starting from 
the lower left corner as before. For a given fixed 777, and large odd q, that is odd q 
which satisfy the condition q > [^^/m + l/2\ + 1, we estimate the elements on the 
corresponding diagonal of X„i_i as follows in view of (|6.21|l 



(6.56) |1 + iiq)\ < 1 + \i{q)\ < 1 + ^ J, - = Ci(7T7) < 2. 



7 y/^ + 1/2 

But for the elements in Xm-i on the diagonals corresponding to small odd q, that 
is for 3 < Q < \-\/ 777 + 1/2| , we have , < — , and for such q we find from 
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Thus for TO > 38 and for q = 3, 5, • • • , 4to — 5 we have 

\l + i{q)\<C{m) 

where C{m) — max(Ci(TO), C2(to)). Now it is elementary to show that for to > 58, 
all the terms on the RHS of (|6.55(l decrease monotonically with to and direct 
computation shows that C2(58) < 1.997 < 2. Together with (|6.56(l this implies 
C(to) < 2 for all to > 58 and hence detT^ 7^ by Proposition 16.11 for to in this 
range. 

But direct evaluation shows that for 38 < to < 57 

V2 to(2to)! / 
and so, again using H6.56(l . we have 

C(to) f I - 4^ • 22™-2) < 1, 38 < TO < 57. 
^ ^2 to(2to)! / ' - - 

Hence det T„i-i ^ also for 38 < to < 57. This completes the proof of Theorem l2.6l 

Appendix 

Our goal here is to prove the two (crude) estimates, 16.44|1 and (|6.51() needed in 
the text. For s G [0, +00), set 



L{s) = s 



dv 



where Rq = 0.855. Note that R{m) > Rq for all jn > 38. Set also 



g — s V V dv 



e 



and 

H{s) = 2(sLi(s) - l). 
Proposition A.l. For s € [0, +00) 

< i(s) < 6. 

Proposition A. 2. We have 



r+oo 

/ \H{s)\ds<2. 
Jo 



The function H{s) changes sign. 

The functions L{s) and H{s) are plotted in Figs. |3| and 01 respectively, using 
Maple. Fig.O "proves" Proposition lA.il Integrating by parts we will prove below 
that for all s > 

(A.2) iy(,)< 1 + ^ + 2y2.4e-V2 

s^ s"^ 

which together with the fact ( "evident" from Fig. 01 and proved below) that H 
is positive on [6, +00) imphes \H{s)\ds < 0.233. Again, using Maple, direct 
numerical integration shows that \H{s)\ds < 2.407. This "proves" Proposition 
I A. 21 with 2.64 on the right. Given the importance of the numerical values of the 
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Figure 3. Graph of L(s) and the line at height 6 




Figure 4. Graph of H{s) 

bounds lt).44|l and Ht).51|l in our Proof of Theorem 12.61 we now give a rigorous 
error analysis (in which we use the computer to the extent of computing the four 
arithmetic operations and the exponential function only) for these bounds. 
In the error analysis both for L and H, the integral 

INT{x)= f e-^''(l- 
Jo 

plays a basic role. Note that for any < a < 1, 

(A.3) / e-^''(l - w)'i/2^v < e"''^2\/r^. 
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Integrating by parts we find 



(A.4) 



2x2 
105 



16a; 



e 



-ay 






-ax ^-j^ 




-3/2j 


-ax ^2 




-5/2j 


-ax ^-j^ 










-9/2 d„ 



1 1 3 15 , , 7/n 

a; 2a;2 4a;3 8x^ 
A.l. Proof of Proposition rOl Clearly L{s) > for s > 0. Note that 



1 1 

e-^^(l-w)-i/2^^> _ 



(A.5) 

X 

Expanding the square of t; = 1 — (1 — w) and integrating by parts we find 



(A.6) 



/' e-^^v'{l ~ v)-'/^ dv^^--^+(l + - + ^) INT{x). 
Jq X 2x^ \ X 4x^/ 



Substututing here the estimates ljA.3|) . ljA.4|) for a = 1/2 and using ljA.5|) we find 
(after simplifications) an upper estimate on (nonnegative) L{s) which implies that 
i(s) < 6 for s > 25. 

Next, making the change of variable t = ^/^ — v in both integrals in L{s) 

L{s) < s^io-f < e^(i-«o) < 6 

/o e-(*^-i) dt - 

for < s < 2. 

Finally consider s G [2,25]. Using (|A.5|I 



1-e-^ Jo y/T^ 

Note that the first factor is < 1.157. Hence it suffices to prove 

L2{s) < 5.185, 2 < s < 25. 

We do this as follows. Assume for a moment that the following (crude) a priori 
bound holds 

(A.7) |L2(s)|<75, 2<s<25. 

Consider for some the mesh Sj = 2+jh, j = 0, 1, • • • , where h — {25~2)/Ne- 
By (HjI) 

(A.8) \L2{s)\<\L2{s,)\+75h, se[sj_i,Sj], j = 1, 2, • • • , TV^. 

The value of L2 at a point is estimated as follows. L2 is s"^ times (jA.6|l with x = Rqs. 



Making the change of variable t = -y/1 — v in INT(x) we obtain for any Ni 
(A.9) INT{Ros) = 2 / e^'"'^*'-^^ '^^ ^ wY. ^^"^ 



Ni 

osiik/Nif-l) 



k=l 
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where the integral is estimated by the right rectangle sum using the monotonia 
growth of e^"**'* forte [0,1]. Now we can just increase A^e, -^i until maxj=i_2.-- , at, 
of the RHS in (|A.8|I (where l|A.9|) has been substituted) becomes < 5.185. It turns 
out that — 8, 000 and Ni ~ 12, 000 suffice and for these values the maximum is 
< 5.162. We must also check that the error does not accumulate too much when we 
sum 12,000 numbers, the precision of each being 10 digits. In the worst scenario, 
the mantissa after each addition accumulates an error of 0.5 in the 10th digit. After 
12,000 additions the error is 6 in the 7th digit. This means that our result has at 
least 5 correct digits. And in the above explanation we have kept only 4 digits in 
arriving at the number 5.162. 

It remains to prove the bound ljA.7|) . Differentiating 6*2(5) we find 

\G',{s)\ < max {s.^ e-«-%^ Ros' e-^-^^ -^}. 



The first integral is estimated using ljA.6|) , (|A.4p . For the second integral: start by 
rewriting the expression in terms of INT{x) (as we have done in (jA.6|l ) and then 
use (|A.4I) . These estimates imply (|A.7|I for s > 2. 

A. 2. Proof of Proposition lA.2l Writing v = 1 — (1 — t;) and integrating by parts 
as above we can rewrite H{s) in terms of INT{s'^) 

(A.IO) i7(,)^ 2.2-1 -—I— = 2.2-1 1 



Combining (|A.10|I with (|A.4|I . the estimate (|A.2p follows. We need also a lower 

estimate on H{s). As in (jA.4|l . for any a E (0, 1) 

(A.11) 



^ dv ^ r dv 1 1 3 15 

XV — ^ I -XV — _ 1 1 1 

~ Jo X 2x2 4x3 82.4 

Together with ljA.10|) . this implies 

> NUMER^jx) 
' - DENOMaix) 

where DENOMaix) is the RHS in (jA.lip and NUMERa{x) is a certain, ex- 
plicit function. It is convenient to take a — 4/5. Then it is elementary to check 
DENOMaix) > for x > 9. We find 

^NUMERM > - ^ - 2v^.=«-'= (1 + I + 5 + ^) 
- 216 17496 

where we have discarded ^ and evaluated the monotonically decreasing terms 
and the term in the paranthesis at x = 9. It is now elementary to check that the 
RHS in (f02)l is > for x > 9. 

We have proved that H{s) > for s > 3. From this and (|A.2|1 it follows that 

/>oo 

(A.13) / \H{s)\ds < 0.233. 



15 



6 
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On the interval [3, 6] where H does not change sign and so \H\ is smooth we use 
the trapezoidal integration scheme with a rigorous error estimation. However on 
[0,3] (where H may and indeed does change sign) we use the following argument. 
Fix Ne and Ni and consider the mesh Sj = jTi, j = 0, 1, • • • , iVg where /i = (3 — 
0)/Ne. Using the left and right rectangles as in ljA.9|l we find that 

HJ<H{s)<H+, j=0,l,---iVe 

where for any Ni 

HT EE 2s' ^ 



3 --"3 ^^{klN{f-l) 



-^"3 + ^ 1 s^iik/N^y-l)- 

This gives 

\H{s)\<ma.^{\Hr\,\H+\}, se[s„Sj+,], j = 0,1, ■ ■ ■ N, - I 

and 

/ ms)\ds<— ^ me.^{\H-\,\H+\} 

for any choice of N.^^,,. Increasing these numbers we improve the estimate. It turns 
out that for ^ N, ^ 3, 000 we get \H{s)\ds < 2.242 (Maple numerical inte- 
gration suggests the estimate 2.200). Again we only use the computer to evaluate 
arithmetic operations and exponents. Also the precision is 10 digits so that af- 
ter 9 • 10^ summations a worst possible error of 0.5 in the 10th digit at each step 
becomes an error of 4.5 in the 4th digit. Thus we write 

(A.14) [ \H {s)\ ds <2.2A7. 

Jo 

Finally we evaluate the integral over [3, 6] (where H is positive) via the trape- 
zoidal scheme to obtain 

r H{s) ds^hJ2 Hisj) + -^^^) + ^(^) h + ERR, + ERR, 
\ERR,\< max |i/"(s)| 

\ERRi\ < (maximal error in computation of H{sj) over all steps) x (6 — 3) 

where Sj = 3 + jh, j — 0,1, - ■ ■ , N,, h ^ {6 — 3)/iVe. At each Sj we will compute 
H{sj) also by the trapezoidal scheme with an explicit error estimate. We demand 
\ERRe\ < 0.05 and \ERR,\ < 0.05. Let us first find the required iVe- Assume for 
a moment a (crude) estimate 

(A.15) |ff"(s)| < 245" -h 205^-1- 4, s > 0. 

Then 

2 

(24 • 6* -I- 2 • 6^ + 4) < 0.05 implies iVe = 1,197. 
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Now consider Ni. For any s £ [3, 6] we compute Fs{u)du where Fs{u) = ^" ^-^ 
(see (|A.10|l ) by the trapezoidal scheme. Note 

52 



max 

ne[o,i] 



+ 2s\ 



Now if Fs {u)du = ^ {s)+erri (s) , then the error in the computation of 1/ Fs {u)dv 
(which enters H{s)) is 



< |em(s)| 



/o ™ " 1/0' I -IE!-)!' 

From e^^ ^'^'^ du ~ ^ e~^^* -?== > — we ffet a lower bound 

min / e"'("'~i)dw > 0.01388. 

ue[3,6] Jo 

When computing the sums below we also store the minimal value of '^{sj) 

over j — 0, 1, • • • , iVe and verify at the end that it is also > 0.01388. Thus we 
demand 

\ERR,\<{6-3) _ max ^''''''^^'^^^ 



j=o,i, -,w. (0.01388)2 



1 f I 



< — max ( — ^— ) 4(4 + s2/2) < 0.05 

- 12i=o,i, -,w. \Ni{sj) J ^ ^ ' (0.01388)2 - 

which holds if we choose 

iV,(s,)> 323^4 + ^2/2, j = 0,l,--- ,A^e 

where Sj = 3 + j . A computer evaluation (again using only elementary opera- 
tions) then shows that /^^ H{s)ds < 0.208 with our guaranteed precision ±0.1. Also 
we keep 10 digits and the total number of summands is < 1, 197- (323-\/6'* + 62/2 + 
1) < 1.5 • 10^. Hence if the worst error of 0.5 in the 10th digit is made at each step, 
the accumulated error does not exceed 7.5 in the 4th digit (in the mantissa of the 
number 2.08- lO^^). Thus 

(A.16) ^ ff(s) ds< 0.208 + 0.1 + 0.00075 = 0.30875 < 0.309. 

(Maple numerical integration gives ^^(s) ds < 0.208.) 
Collecting (|A.13p . (|A.14p . (|A.16|) we conclude 

/•OO 

/ \H{s)\ ds < 0.233 + 2.247 + 0.309 = 2.789 < 2.8 
Jo 

which completes the proof up to (|A.15|) . To prove l|A.15|l denote 

K{x)^-2 + 2x-{^, /^'^ /'e— ^iL, Z = 0,l,2,.... 
J 1 J Jo V^-v 

Then 

(A.17) His)=K{s^), H'{s)^2s-K'is^). 
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We find 



K'{x) = 2-^ + 2x- 




If we single out the positive and the negative terms and use the fact that each 
fraction is bounded from above by 1, we obtain 

-2x < K'{x) < 2 + 2x, a; > 0. 

and from IIA.171) 




|iJ'(s)| < 4/ + 4s, s>0. 
After computing H" {s) in terms of K'(x), K" [x) in the same way, we arrive at 
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